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Abstract. Three-dimensional hydrodynamical, Newto- 
nian calculations of the coalescence of equal-mass binary 
neutron stars are performed with the "Piecewise Parabolic 
Method" . The properties of neutron star matter are de- 
scribed by the equation of state of Lattimer & Swesty 
(1991) which allows us to include the emission of neutri- 
nos and to evaluate our models for the ^^-annihilation in 
the vicinity of the merging stars. When the stars have 
merged into one rapidly spinning massive body, a hot 
toroidal cloud of gas with a mass of about 0.1-0.2 M© 
forms around the wobbling and pulsating central ~ 3 Mq 
object. At that time the total neutrino luminosity climbs 
to a maximum value of 1-1.5 • 10 53 erg/s of which 90-95% 
originate from the toroidal gas cloud surrounding the very 
dense core. The mean energies of u e , D e , and heavy- lepton 
neutrinos v x are around 12 MeV, 20 MeV, and 27 MeV, 
respectively. The characteristics of the neutrino emission 
are very similar to the emission from type-II supernovae, 
except for the D e luminosity from the merged neutron stars 
which is a factor 3-6 higher than the luminosities of the 
other neutrino species. 

When the neutrino luminosities are highest, vv- 
annihilation deposits about 0.2-0.3% of the emitted neu- 
trino energy in the immediate neighborhood of the merger, 
and the maximum integral energy deposition rate is 3- 
4 ■ 10 50 erg/s. Since the 3 Mq core of the merged object 
will most likely collapse into a black hole within millisec- 
onds, the energy that can be pumped into a pair-photon 
fireball is insufficient by a factor of about 1000 to explain 
7-ray bursts at cosmological distances with an energy of 
the order of 10 51 /(47r) erg/steradian. Analytical estimates 
show that the additional energy provided by the annihi- 
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lation of vv pairs emitted from a possible accretion torus 
of ~ 0.1 Mq around the central black hole is still more 
than a factor of 10 too small, unless focussing of the 
fireball into a jet-like expansion plays an important role. 
A few 10 -4 Mq of very neutron-rich, low-entropy matter 
may be dynamically ejected shortly after the neutron stars 
have merged, and another 10~ 4 up to a few 10~ 2 Mq of 
strongly neutronized, high-entropy material could be car- 
ried away from the accretion torus in a neutrino-driven 
wind. The contamination with this baryonic material is a 
severe threat to a relativistic fireball. Aspects of a possible 
r-processing in these ejecta are discussed. 

Key words: gamma rays: bursts - elementary particles: 
neutrinos nuclear reactions, nucleosynthesis, abundances 
- stars: neutron - binaries: close - hydrodynamics 



1. Introduction 

During the first two and a half years of operation, the 
Burst and Transient Source Experiment (BATSE) on the 
Compton Gamma Ray Observatory (GRO) has observed 
1122 cosmic gamma-ray bursts (Meegan et al. 1995a). The 
angular distribution of the bursts on the sky is amazingly 
well compatible with isotropy (e.g., Tegmark et al. 1995). 
There are recent claims that a tenuous indication of an as- 
sociation with host galaxies is present in the burst sample 
evaluated by Larson et al. (1996) and that a correlation 
with Abell clusters has been found on a 95% confidence 
level in the BATSE 3B catalog (Kolatt & Piran 1996). If 
confirmed, these would be the first hints of a correlation 
of gamma-ray bursts with any other astronomical popu- 
lation. Number counts show a deficiency of faint bursts 
relative to the —3/2 power law expected for the bright- 
ness distribution of a homogeneous spatial distribution of 
standard candle sources (Fenimore 1996). This paucity of 
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weak bursts could be caused by a truncation of the dis- 
tance to very far burst sources or might be an effect due to 
the expansion of the universe or could be associated with 
an evolution of the spatial density of bursts (Bloom et 
al. 1996). Both observational facts, isotropy and bend in 
the brightness distribution, would be naturally explained 
if the gamma-ray bursters were situated at cosmological 
distances (e.g., Paczyhski 1995, Hartmann et al. 1996, Ko- 
latt & Piran 1996). Nevertheless, the assumption that they 
might populate an extended Galactic halo cannot be ruled 
out yet (Lamb 1995), and Galactic halo models have been 
constructed (Podsiadlowski et al. 1995) which are able to 
fulfill the stringent limits set by the isotropy of the de- 
tected bursts and their inhomogeneous spatial distribu- 
tion. 

The distribution of measured gamma-ray burst dura- 
tions exhibits a bimodal structure with peaks at about 
0.5 s and about 30 s (Meegan et al. 1995a). The bursts can 
be as short as ~ 1 ms but can also last for several 100 s 
with variabilities and fluctuations on a millisecond time 
scale (Fishman et al. 1994). The extremely short sub-ms 
rise times of the gamma-ray luminosities suggest that the 
energy sources for the bursts must be connected with very 
compact astrophysical objects which have a typical size of 
the order of 100 km. This favors neutron stars or black 
holes as most likely candidates for the enigmatic origin of 
the cosmic gamma-ray bursts. 

Despite of more than 25 years of gamma-ray burst ob- 
servations, there is neither a convincing identification of 
counterparts in any other energy range of the electromag- 
netic spectrum (Greiner 1995a, 1995b; Vrba 1996), nor has 
a generally accepted, satisfactory theoretical model been 
developed yet (Nemiroff 1994a, Hartmann & Woosley 
1995, Woosley 1996). About 120 gamma-ray burst mod- 
els have been published in the refereed literature until 
1992 (Nemiroff 1994a), until 1994 there were 135 (Nemiroff 
1994b), and maybe another one or two dozens have been 
added since. 

Cosmological explanations have become increasingly 
popular in the more recent publications, a fair fraction 
of which suggests collisions of two neutron stars or merg- 
ers of binaries consisting of cither two neutron stars (NS- 
NS) or a black hole and a neutron star (BH-NS) as pos- 
sible sources of the bursts (e.g., Paczyflski 1986; Good- 
man 1986; Eichler et al. 1989; Piran 1990; Paczynski 1991; 
Narayan et al. 1991, 1992; Piran et al. 1992; Meszaros 
& Rees 1992a,b, 1993; Woosley 1993a; Mochkovitch et 
al. 1993, 1995a; Hcrnanz et al. 1994, Katz & Canel 1995a). 
One of the reasons for the attractivity of these scenarios 
is the desired compactness of the objects, another rea- 
son is the knowledge that these events should happen 
and should release large amounts of energy (Dermer & 
Weiler 1995). The frequency of NS-NS and BH-NS merg- 
ers was estimated to be between 10~ 6 and 10~ 4 per year 
per galaxy (Narayan et al. 1991, Phinney 1991, Tutukov et 
al. 1992, Tutukov & Yungelson 1993, Lipunov et al. 1995) 



and is therefore sufficient to explain the observed burst 
rate which requires an event rate of about 10~ 6 yr^ 1 per 
galaxy (Narayan et al. 1992). These rates per galaxy are 
so low that burst repetition in the same region of the sky 
is practically excluded which is in agreement with the ob- 
servations (Lamb 1996, Meegan et al. 1995b, Brainerd et 
al. 1995, Efron & Petrosian 1995). If the merger rate is 
near the high end of the estimated range, some beam- 
ing of the gamma-ray emission might be involved, or the 
majority of the bursts has to be very dim and escapes de- 
tection. Beaming would also lower the energy that must 
be converted into gamma rays at the source in order to 
cause the observed fluences. 

Since the detected gamma-ray bursts appear to be 
isotropically distributed and do not visibly trace the large- 
scale structure of luminous matter in the universe, in par- 
ticular, are not concentrated towards the supergalactic 
plane like nearby galaxies, constraints on the distance 
scale to cosmological bursts can be placed. From the 
BATSE 3B catalog Quashnock (1996) infers that the co- 
moving distance to the "edge" of the burst distribution is 
greater than 630 hr x Mpc and the nearest bursts are far- 
ther than AOh^ 1 Mpc (at the 95% confidence level), the 
median distance to the nearest burst being 170 h^ 1 Mpc 
(h is the Hubble constant in units of 100 km/s/Mpc). 
From the absence of anisotropics in supergalactic coordi- 
nates, Hartmann et al. (1996) conclude that the minimum 
sampling distance is 200 ft. -1 Mpc, and Kolatt & Piran 
(1996) find for their accurate position sub-sample mem- 
bers locations within 600 h^ 1 Mpc. In case of isotropic 
emission, standard candle non-evolving burst sources at 
these cosmological distances must release 7-ray energies 
of the order of (3. ..4) • lO 51 ^ 2 erg (Woods & Loeb 1994, 
Quashnock 1996). 

This energy is about 0.1-0.2% of the rest-mass energy 
of one solar mass or roughly 1% of the gravitational bind- 
ing energy set free when two 1.5 Mq neutron stars merge. 
A large part of the energy released during the merging, i.e., 
up to more than 10% of M Q c 2 , is carried away by gravi- 
tational waves, the exact value depending on the nuclear 
equation of state and thus on the compactness of the neu- 
tron stars and of the merged object (see Ruffert et al. 1996 
and references therein) . A similar amount of energy could 
be radiated away in neutrinos which are abundantly pro- 
duced when the matter of the coalescing and merging stars 
is heated up to very high temperatures by tidal forces, fric- 
tion, and viscous dissipation of kinetic energy in shocks 
(Eichler et al. 1989, Narayan et al. 1992, Harding 1994). 
The duration of the neutrino emission, the neutrino lu- 
minosity, and the total energy radiated in neutrinos will 
be determined by the structure and dynamical evolution 
of the merger, by the thermodynamical conditions in the 
merging stars, and by the lifetime of the merged object 
before it collapses into a black hole or before the surround- 
ing material is swallowed by the black hole. Gravitational 
waves as well as neutrinos from these very distant sources 
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cannot be measured with current experiments, but detec- 
tors are planned or built for gravitational waves (LIGO, 
VIRGO, GEO600; see, e.g., Thorne 1992) and are envi- 
sioned for neutrinos (a kilometer-scale neutrino telescope 
with 1 km 3 of instrumented volume; see Weiler et al. 1994, 
Halzen & Jaczko 1996, Halzen 1996), which might be lucky 
to catch the death throes of massive binaries in the not 
too distant future. 

Due to the very high opacity of neutron star matter, 
high energy photons cannot be emitted by the merging 
object directly, unless the very outer layers are heated to 
sufficiently high temperatures. Even if they were radiat- 
ing with luminosities substantially above the Eddington 
limit (Duncan et al. 1986), NS-NS or NS-BH mergers at 
cosmological distances would still be far too faint to be 
visible from Earth. However, if only a tiny fraction (less 
than 1%) of the potentially emitted neutrinos and antineu- 
trinos annihilate in the vicinity of the merger (Goodman 
et al. 1987, Cooperstein et al. 1987) and create a fireball 
of electron-positron pairs and photons (Cavallo & Rccs 
1978), an intense outburst of gamma-radiation can be pro- 
duced with an overall power output exceeding the Edding- 
ton luminosity by up to 15 orders of magnitude and en- 
ergetic enough to account for a cosmological gamma-ray 
burst (Eichler et al. 1989). 

Relativistic expansion of the pair-photon plasma and 
the final escape of high energy gamma-radiation with the 
observed non-thermal spectrum (e.g., Band 1993) from an 
optically thin fireball can only occur if the baryon load of 
the radiation-pressure ejected shells is sufficiently small, 
i.e., the contamination with baryons must be below a cer- 
tain limit (Goodman 1986, Paczynski 1986). If the ad- 
mixture of baryons is too large, highly relativistic outflow 
velocities will be impeded through conversion of radiation 
energy to kinetic energy and the flow will remain optically 
thick, leading to a degradation of the photons to the UV 
range (Paczynski 1990, Shemi & Piran 1990). In order to 
suppress pair production 7 + 7 ^ e+ + e~ and to make the 
fireball optically thin to its own photons, the Lorentz fac- 
tor r of the emitting region must obey T 100 (Fenimore 
et al. 1993, Meszaros et al. 1993). The bulk Lorentz factor 
is related with the ratio of total energy to (baryon) rest 
mass energy: T = E 1 / (£Mc 2 ) where £ is the fraction of the 
total energy in the fireball that ends up as observed energy 
in 7-rays, Ej. Therefore a lower limit for T places an upper 
bound on the amount of baryonic mass M released in the 
explosive event or being present as ambient gas near the 
burst source, M <, 0.6 • lO^^H^/iO 51 erg) M Q . Based 
on the empirical data from the first BATSE catalog (169 
bursts), Woods & Locb (1994) deduce a mean minimum 
Lorentz factor T of about 500, corresponding to an aver- 
age maximum baryon load of 10 -5 £ -1 M© for gamma-ray 
burst events at cosmological distances. Recent introduc- 
tory reviews and summaries of the properties of fireballs in 
the context of cosmological models of gamma-ray bursts 



were given by Dermer & Weiler (1995), Meszaros (1995), 
and Piran (1995, 1996). 

The limits for the baryon load in the fireball set by ob- 
servational requirements are very stringent and the bary- 
onic pollution of the pair-photon plasma is a major con- 
cern for gamma-ray burst models based on mergers of mas- 
sive binary stars. In order to obtain ^-annihilation in a 
baryon-poor region around the merger, anisotropies of the 
merger geometry are considered to be essential (Narayan 
et al. 1992; Meszaros & Rees 1992a,b; Woosley 1993a; 
Mochkovitch 1993, 1995a; Hcrnanz 1994). The hope is that 
centrifugal forces can keep a region near the rotation axis 
of the merging NS-NS or NS-BH system relatively baryon 
free and that the neutrino-driven mass loss from an accre- 
tion disk or torus formed after the merging of the binary 
will leave a "clean" funnel along the axis of symmetry 
where ^-annihilation can create collimated e + e~ jets ex- 
panding relativistically in both axis directions. Also, gen- 
eral relativistic bending of the v and v trajectories and vv- 
annihilation inside the innermost stable orbit for massive 
particles around the accreting black hole have been sug- 
gested to be potentially helpful. Secondary processes could 
lead to the reconversion of baryonic kinetic energy, e.g. by 
external shock interactions of the expanding fireball in the 
ambient interstellar medium or pre-ejected gas (Rees & 
Meszaros 1992; Meszaros & Rees 1992a, 1993; Sari & Pi- 
ran 1995) or by internal shock interactions of shells having 
different speeds within the expanding fireball (Narayan et 
al. 1992, Paczynski & Xu 1994, Rees & Meszaros 1994, 
Mochkovitch et al. 1995b). Even more complex physics 
like the upscattering of interstellar photons of a local ther- 
mal radiation field by collisions with electrons in the ul- 
trarelativistic wind (Shemi 1994) or highly amplified mag- 
netic fields (Meszaros & Rees 1992a, Usov 1992, Smolsky 
& Usov 1996) might be important in determining the tem- 
poral and spectral structure of the observable gamma-ray 
burst. In fact, such secondary processes may be crucial 
to produce the detected very high energy photons with 
energies of up to more than 1 GeV (Hurley et al. 1994, 
Teegarden 1995). 

In this work we shall not deal with the possibly very 
complex and complicated processes that are involved in 
the formation of the finally observable gamma-ray sig- 
nature. Instead, our interest will be concentrated on the 
hydrodynamical modelling of the last stages of the coa- 
lescence of binary neutron stars employing an elaborate 
equation of state for neutron star matter with the aim to 
compute the neutrino emission from the merger. The re- 
sults of our models will allow us to calculate the efficiency 
of ^-annihilation and the energy deposition in the vicin- 
ity of the merging stars. After dozens of papers that sug- 
gest and refer to the annihilation of vv pairs from merg- 
ing compact binaries as the energy source of cosmologi- 
cal gamma-ray bursts, we shall try to put this hypothesis 
to a quantitative test. We shall investigate the questions 
whether the neutrino emission from the tidally heated neu- 
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tron stars just prior to merging (Meszaros & Rees 1992b), 
during the dynamical phase of the merging or collision 
(Narayan et al. 1992, Meszaros & Rees 1992b, Dermer & 
Weiler 1995, Katz & Canel 1995a) or after the merging 
when a hot accretion disk or torus has possibly formed 
around a central black hole (Woosley 1993a, Mochkovitch 
et al. 1993), is sufficiently luminous and lasting to yield 
the energy required for gamma-ray bursts at cosmological 
distances. Also, our simulations will provide information 
about how much mass might remain in an accretion disk 
or torus and about the thermodynamcial conditions in this 
disk matter. These aspects might have interesting implica- 
tions for the possible contributions of NS-NS and NS-BH 
mergers to the nucleosynthesis of heavy elements. 

The paper is organized as follows. In Sect. ^ the com- 
putational method and initial conditions for our simula- 
tions are described and the hydrodynamical evolution of 
the merger is shortly summarized from the results given 
in detail by Ruffert et al. 1996 (Paper I). In Sect. ||the re- 
sults for the neutrino emission and thermodynamical evo- 
lution of the merger will be presented. Section |] deals with 
the neutrino- antincutrino annihilation and contains infor- 



mation about the numerical evaluation (Sect. 4.1), about 
the numerical results (Sect. |4~2] ), and about an analyti- 
cal model which was developed to estimate the neutrino 
emission and annihilation for an accretion torus around 
a black hole (Sect. 4.3). In Sect. || the results will be dis- 
cussed concerning their implications for heavy element nu- 
cleosynthesis and for gamma-ray bursts, and Sect. ^| con- 
tains a summary and conclusions. 

2. Computational procedure, initial conditions, 
and hydrodynamical evolution 

In this section we summarize the numerical methods and 
the treatment of the input physics used for the presented 
simulations. In addition, we specify the initial conditions 
by which our different models are distinguished. Also, the 
results for the dynamical evolution as described in detail 
in Paper I are shortly reviewed. 

2.1. Numerical treatment 

The hydrodynamical simulations were done with a code 
based on the Piecewise Parabolic Method (PPM) devel- 
oped by Colella k Woodward (1984). The code is basi- 
cally Newtonian, but contains the terms necessary to de- 
scribe gravitational wave emission and the corresponding 
back-reaction on the hydrodynamical flow (Blanchet et 
al. 1990). The modifications that follow from the grav- 
itational potential are implemented as source terms in 
the PPM algorithm. The necessary spatial derivatives are 
evaluated as standard centered differences on the grid. 

In order to describe the thermodynamics of the neu- 
tron star matter, we use the equation of state (EOS) of 
Lattimer & Swesty (1991) in a tabular form. The inver- 



sion for the temperature is done with a bisection iteration. 
Energy loss and changes of the electron abundance due to 
the emission of neutrinos and antineutrinos are taken into 
account by an elaborate "neutrino leakage scheme". The 
energy source terms contain the production of all types 
of neutrino pairs by thermal processes and of electron 
neutrinos and antineutrinos also by lepton captures onto 
baryons. The latter reactions act as sources or sinks of 
lepton number, too, and are included as source term in a 
continuity equation for the electron lepton number. When 
the neutron star matter is optically thin to neutrinos, the 
neutrino source terms are directly calculated from the re- 
action rates, while in case of optically thick conditions 
lepton number and energy are released on the correspond- 
ing neutrino diffusion time scales. The transition between 
both regimes is done by a smooth interpolation. Matter 
is rendered optically thick to neutrinos due to the main 
opacity producing reactions which are neutrino-nucleon 
scattering and absorption of neutrinos onto baryons. 

More detailed information about the employed numeri- 
cal procedures can be found in Paper I, in particular about 
the implementation of the gravitational wave radiation 
and back-reaction terms and the treatment of the neutrino 
lepton number and energy loss terms in the hydrodynam- 
ical code. 

2.2. Initial conditions 

We start our simulations with two identical Newtonian 
neutron stars with a baryonic mass of about 1.63 Mq and 
a radius of 15 km which are placed at a center-to-center 
distance of 42 km on a grid of 82 km side length. With the 
employed EOS of Lattimer & Swesty (1991), this baryonic 
mass corresponds to a (general relativistic) gravitational 
mass of approximately 1.5 Mq for a cool star with a ra- 
dius as obtained from the general relativistic stellar struc- 
ture equations of 11.2 km. Having a compressibility mod- 
ulus of bulk nuclear matter of K = 180 MeV (which is 
the "softest" of the three available cases) , the Lattimer & 
Swesty EOS may overestimate the stiffness of supranuclear 
matter, in particular, since in its present form it neglects 
the possible occurrence of new hadronic states besides the 
neutron and proton at very high densities. For a softer 
supranuclear EOS, neutron stars would become more com- 
pact and their gravitational binding energy larger so that 
a baryonic mass of 1.63 Mq would more likely correspond 
to a gravitational mass between 1.4 and 1.45 Mq. 

The distributions of density p and electron fraction Y e 
are taken from a one-dimensional model of a cold, delep- 
tonized neutron star in hydrostatic equilibrium. For nu- 
merical reasons the surroundings of the neutron stars can- 
not be assumed to be evacuated. The density of the ambi- 
ent medium was set to 10 9 g/cm 3 , more than five orders of 
magnitude smaller than the central densities of the stars. 
The internal energy density and electron fraction of this 
gas were taken to be equal to the values in the neutron 
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stars at a density of 10 9 g/cm 3 . In order to ensure suffi- 
ciently good numerical resolution, we artificially softened 
the extremely steep density decline towards the neutron 
star surfaces by not allowing for a density change of more 
than two orders of magnitude from zone to zone. From 
this prescription we obtain a thickness of the neutron star 
surface layers of about 3 zones. 



Table 1. Parameters and some computed quantities for all 
models. N is the number of grid zones per dimension in the or- 
bital plane, S defines the direction of the spins of the neutron 
stars relative to the direction of the orbital angular momen- 
tum, and fceTex gives the maximum temperature (in energy 
units) reached on the grid during the simulation of a model. 
Z/„ c is the electron neutrino luminosity after approaching a sat- 
uration level at about 8 ms, L^ c is the corresponding electron 
antineutrino luminosity, and L Vx the luminosity of each indi- 
vidual species of u^, P M , u T , v T . L v gives the total neutrino 
luminosity after a quasi-stationary state has been reached at 
t Ji 6-8 ms and E vii denotes the integral rate of energy depo- 
sition by neutrino-antineutrino annihilation at that time 



of the numerical resolution of the computations with 64 3 
zones. 

The rotational state of the neutron stars is determined 
by the action of viscosity. If the dynamic viscosity of neu- 
tron star matter were large enough, tidal forces could lead 
to tidal locking of the two stars and thus spin-up dur- 
ing inspiral. Kochanek (1992), Bildsten & Cutler (1992), 
and Lai (1994), however, showed that microscopic shear 
and bulk viscosities are probably orders of magnitude too 
small to achieve corotation. Moreover, they argued that 
for the same reason it is extremely unlikely that the stars 
are heated up to more than about 10 8 K by tidal inter- 
action prior to merging. In this sense, models A64 and 
A128 can be considered as reference case for two non- 
corotating neutron stars, while model B64 represents the 
case of rigid-body like rotation of the two stars, and model 
C64 the case where the spin directions of both stars were 
inverted. 

Since in the case of degenerate matter the tempera- 
ture is extremely sensitive to small variations of the in- 
ternal energy, e.g. caused by small numerical errors, we 
did not start our simulations with cold (T = 0) or "cool" 
(T ^ 10 s K) neutron stars as suggested by the investiga- 



Model 


N 


S 










L v 








McV 


10 53 crg/s 


10 53 org/s 


10 53 erg/s 


10 53 crg, 


A64 


64 





40 


0.13 


0.58 


0.11 


1.15 


B64 


64 


+1 


30 


0.18 


0.55 


0.09 


1.09 


C64 


64 


-1 


>50 


0.15 


0.67 


0.10 


1.22 


A128 


128 





39 


0.16 


0.43 


0.06 


0.83 



/s 10 



The orbital velocities of the coalescing neutron stars 
were prescribed according to the motion of point masses 
spiralling towards each other due to the emission of grav- 
itational waves. The tangential velocities of the neutron 
star centers were set equal to the Kepler velocities on 
circular orbits and radial velocity components were at- 
tributed as calculated from the quadrupole formula. In 
addition to the orbital angular momentum, spins around 
their respective centers were added to the neutron stars. 
The assumed spins and spin directions were varied be- 
tween the calculated models. Table |l| lists the distinguish- 
ing model parameters, the number N of grid zones per 
spatial dimension (in the orbital plane) and the spin pa- 
rameter S. The neutron stars in models A64 and A128 
did not have any additional spins (S = 0), in model B64 
the neutron star spins were parallel to the orbital angu- 
lar momentum vector (S — +1), in model C64 the spins 
were in the opposite direction (S = —1). In both models 
B64 and C64, the angular velocities of the rigid neutron 
star rotation and of the orbital motion were chosen to be 
equal. Model A128 had the same initial setup as model 
A64 but had twice the number of grid zones per spatial 
dimension and thus served as a check for the sufficiency 



lotions of Kochanek (1992), Bildsten & Cutler (1992), and 
3Q CI ^aj (1994). Instead, we constructed initial temperature 

2 ^distributions inside the neutron stars by assuming ther- 

3 (jrial energy densities of about 3% of the degeneracy energy 
3 ^density for a given density p and electron fraction Y e . The 
— corresponding central temperature was around 7 MeV and 

the average temperature was a few MeV and thus of the or- 
der of the estimates obtained by Meszaros & Rees (1992b) 
for the phase just prior to the merging. Locally, these ini- 
tial temperatures were much smaller than the tempera- 
tures produced by the compression and dissipative heat- 
ing during coalescence (see Table [l] for the maximum tem- 
peratures). However, they are orders of magnitude larger 
than can be achieved by tidal dissipation with plausible 
values for the microscopic viscosity of neutron star mat- 
ter. Even under the most extreme assumptions for viscous 
shear heating, the viscosity of neutron star matter turns 
out to be at least four orders of magnitude too small (see 
Janka & Ruffert 1996). 

Models A64, B64, and C64 were computed on a Cray- 
YMP 4/64 where they needed about 16 MWords of main 
memory and took approximately 40 CPU-hours each. For 
the better resolved model A128 we employed a grid with 
128 x 128 x 64 zones instead of the 64 x 64 x 32 grid of 
models A64, B64, and C64. Note that in the direction or- 
thogonal to the orbital plane only half the number of grid 
zones was used but the spatial resolution was the same 
as in the orbital plane. Model A128 was run on a Cray- 
EL98 4/256 and required about 22 MWords of memory 
and 1700 CPU-hours. 
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2.3. Hydrodynamical evolution 

The hydrodynamical evolution and corresponding gravi- 
tational wave emission were detailed in Paper I. Again, 
we only summarize the most essential aspects here. 

The three-dimensional hydrodynamical simulations 
were started at a center-to-center distance of the two neu- 
tron stars of 2.8 neutron star radii. This was only slightly 
larger than the separation where the configurations be- 
come dynamically unstable which is at a distance of ap- 
proximately 2.6 neutron star radii. Gravitational wave 
emission leads to the decay of the binary orbit, and al- 
ready after about one quarter of a revolution, approx- 
imately 0.6 ms after the start of the computations, the 
neutron star surfaces touch because of the tidal deforma- 
tion and stretching of the stars. 

When the two stars begin to plunge into each other, 
compression and the shearing motion of the touching sur- 
faces cause dissipation of kinetic energy and lead to a 
strong increase of the temperature. From initial values of 
a few MeV, the gas heats up to peak temperatures of sev- 
eral ten MeV. In case of model C64, a temperature of 
nearly 50 MeV is reached about 1 ms after the stars have 
started to merge. After one compact, massive body has 
formed from most of the mass of the neutron stars about 
3 ms later, more than 50 MeV are found in two distinct, 
extremely hot off-center regions. 

At the time of coalescence and shortly afterwards, 
spiral-arm or wing-like extensions are formed from mate- 
rial spun off the outward directed sides of the neutron stars 
by tidal and centrifugal forces. Due to the retained angu- 
lar momentum, the central body performs large-amplitude 
swinging motions and violent oscillations. This wobbling 
of the central body of the merger creates strong pressure 
waves and small shocks that heat the exterior layers and 
lead to the dispersion of the spiral arms into a vertically 
extended "ring" or thick toroidal disk of gas that sur- 
rounds the massive central object. While the mass of the 
compact body is larger than 3 M© and its mean density 
above 10 14 g/cm 3 , the surrounding cloud contains only 
0.1-0.2 M & and is more dilute with an average density 
of about 10 12 g/cm 3 . 

The central body is too massive to be stabilized by gas 
pressure for essentially all currently discussed equations 
of state of nuclear and supranuclear matter. Moreover, 
it can be argued (see Paper I) that its angular momen- 
tum is not large enough to provide rotational support. 
Therefore, we expect that in a fully general relativistic 
simulation, the central object will collapse into a black 
hole on a time scale of only a few milliseconds after the 
neutron stars have merged. Dependent on the total angu- 
lar momentum corresponding to the initial setup of the 
neutron star spins, some of the material in the outer re- 
gions of the disk obtains enough angular momentum (e.g., 
by momentum transfer via pressure waves) to flow across 
the boundary of the computational grid at a distance of 



about 40 km from the center. In model B64 which has 
the largest total angular momentum because of the as- 
sumed solid-body type rotation, 0.1-0.15 M Q are able to 
leave the grid between 2 ms and 4 ms after the start of 
the simulation. However, at most a few times 10 -4 Mq of 
this material have a total energy that might allow them 
to escape from the gravitational field of the merger. In 
model C64 the anti-spin setup leads to violent oscillations 
of the merged body which create a very extended and peri- 
odically contracting and reexpanding disk. In this model 
nearly 0.1 Mq are lost across the outer grid boundaries 
even during the later stages of the simulation (t *S 5 ms). 

While by far the major fraction of the merger mass 
(^ 3 Mq) will be swallowed up by the forming black 
hole almost immediately on a dynamical time scale, some 
matter with sufficiently high specific angular momentum 
might be able to remain in a disk around the slowly rotat- 
ing black hole. Estimates show (Paper I) that with the 
typical rotation velocities obtained in our models, this 
matter must orbit around the central body at radii of at 
least between about 40 km and 100 km. This means that 
only matter that has been able to leave the computational 
grid used in our simulations has a chance to end up in a 
toroidal disk around the black hole. From these arguments 
we conclude that such a possible disk might contain a mass 
of at most 0.1-0.15 Mq. 

3. Neutrino emission and thermodynamical evolu- 
tion 

3.1. Neutrino emission 

The local energy and lepton number losses due to neu- 
trino emission are included via source terms in our code 
as described in the appendix of Paper I. We treated the 
neutrino effects in terms of an elaborate leakage scheme 
that was calibrated by comparison with results from diffu- 
sion calculations in one-dimensional situations. By adding 
up the local source terms over the whole computational 
grid, one obtains the neutrino luminosities of all individual 
neutrino types, the sum of which gives the total neutrino 
luminosity. In the same way number fluxes of electron neu- 
trinos (v e ), electron antineutrinos {y e ), and heavy- lepton 
neutrinos (y^, P^, v T , and P T , which will be referred to as 
v x in the following) can be calculated. The mean energies 
of the emitted neutrinos result from the ratios of neutrino 
luminosities to neutrino number fluxes. 

Figures |l[^ display the results for the neutrino fluxes 
and mean energies of the emitted neutrinos for our mod- 
els A64, B64, C64, and A128. In Fig. |l| one can see that 
the total neutrino luminosities start to increase above 
~ 10 52 erg/s at about 2.5-3.5 ms after the start of the 
simulations. This is about the time when the spiral arm 
structures become dispersed into a spread-out cloud of 
material that surrounds the merger and is heated by the 
interaction with compression waves. At the same time the 
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Fig. 1. Total neutrino luminosities as functions of time for all 
four models 



Fig. 2. Comparison of the cumulative energies emitted in grav- 
itational waves and in neutrinos as functions of time for all four 
models 




Fig. 3. Luminosities of individual neutrino types (v e , v e , and 
the sum of all v x ) as functions of time for the three models 
A64, B64, and C64 



Fig. 4. Average energies of emitted neutrinos v e , v e , and v. 
as functions of time for the three models A64, B64, and C64 



temperatures in the interior of the massive central body 
reach their peak values. Since the neutrino luminosities 
are by far dominated by the contributions from the disk 
emission (see below), the heating of the cloud and torus 
material is reflected in a continuous increase of the neu- 
trino fluxes. Also, the dynamical expansions and contrac- 
tions caused by the oscillations and wobbling of the central 
body impose fluctuations on the light curves. The peri- 
odic expansions and contractions have particularly large 



amplitudes in model C64 because of the anti-spin setup 
of the neutron star rotations which leads to very strong 
internal shearing and turbulent motions after merging as 
well as to higher temperatures than in the other models 
(see Table |). The light curve fluctuations in model C64 
are as large as 25-30% of the average luminosity and pro- 
ceed with a period of 0.5-1 ms which is about the dynam- 
ical time scale of the merger. When the oscillating central 
body enters an expansion phase, very hot matter that is 
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Fig. 5. Energy emission rates (in erg/cm 3 /s) of electron neutrinos (panel a), electron antineutrinos (panel b), the sum of all 
heavy-lepton neutrinos (panel c) , and the total neutrino energy loss rate in the orbital plane of model A64 at the end of the 
simulation (time in the top right corner of the panels). The contours are logarithmically spaced in intervals of 0.5 dex, bold 
contours are labeled with their respective values. The grey shading emphasizes the emission levels, dark grey corresponding to 
the strongest energy loss by neutrino emission 



M. Ruffert et al.: Coalescing neutron stars - a step towards physical models 



9 



Fig. 6. Same as Fig. |5| but for model A128 at time t — 8.80 ms. The higher resolution of this simulation allows more fine 
structure to be visible 
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Fig. 7. Energy loss rates by neutrino emission in two orthogonal planes vertical to the orbital plane for model A128 at time 
t = 8.80 ms. The displayed information is the same as in Fig. 
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Fig. 8. Same as Fig. g but for model B64 at time t = 10.55 ms 

located in a shell around the central core of the merger, is 
swept to larger radii. In course of the expansion the neu- 
trino optical depth decreases and the very hot material 
releases enhanced neutrino fluxes. The neutrino outburst 
is terminated, when the surface-near matter has cooled 
by adiabatic expansion or when re-contraction sets in and 
the neutrino optical depth increases again. At times later 



than about 6-8 ms quasi-stationary values of the fluxes 
are reached which are between about 8 • 10 52 erg/s in case 
of model A128 and about 1.3 • 10 53 erg/s for C64. 

Figure || shows that the energy lost in gravitational 
waves during the merging is about two orders of mag- 
nitude larger than the energy radiated away in neutri- 
nos during the simulated period of approximately 10 ms. 
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Fig. 9. Same as Fig. g but for model C64 at time t — 10.86 ms 

While the gravitational wave luminosity peaks around 
the time when the dynamical instability of the orbit sets 
in and the two neutron stars start to interact dynam- 
ically and fuse into a single object (t between 0.5 ms 
and 1.5 ms), the energy emitted in neutrinos becomes siz- 
able only after the extended toroidal cloud of matter has 
formed around the merged stars. 



The heated torus or "disk" consists of decompressed 
neutron star matter with an initially very low electron 
number fraction Y e between about 0.01 and 0.04. The neu- 
trino emission of the disk is therefore clearly dominated 
by the loss of electron antineutrinos which are primarily 
produced in the process e + + n — > p+v e , because positrons 
are rather abundant in the hot and only moderately de- 
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Fig. 10. Same as Fig. M but for model C64 at time t — 10.86 ms 



generate, neutron-rich matter (see Sect. 3.2). In Fig. gone 
sees that the v e luminosity L Pc is a factor 3-4.5 larger than 
L Ve and between 5.5 • 10 52 erg/s for model B64 and about 
6.7 • 10 52 erg/s for C64 when quasi-stationary conditions 
have been established after about 6 ms from the start of 
the simulations (Table |l|) . The sum of all heavy- lepton 
neutrino fluxes which is four times the individual lumi- 



nosities of v^, Pfj,, v r or 9 T , is around 4 • 10 52 erg/s. The 
better resolved model A128 yields slightly smaller values 
for all luminosities than A64 (Table |l|). 

The mean energies of the emitted neutrinos are dis- 
played in Fig. | for models A64, B64, and C64. The dif- 
ferences between the different models are smaller than in 
case of the neutrino luminosities. This means that the 
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effective temperatures in the neutrinospheric regions are 
very similar in all models. Electron neutrinos are emitted 
with an average energy of 12-13 MeV, electron antineu- 
trinos have 19-21 MeV, and heavy lepton neutrinos leave 
the merger with mean energies of 26-28 MeV during the 
stationary phase. Except for the dominance of relative 
to Z/„ e , both the mean neutrino energies and the neutrino 
luminosities are in good overall agreement with typical 
numbers obtained in stellar core collapse and supernova 
simulations at a stage some time after the prompt v e burst 
has been emitted and after the collapsed stellar core has 
heated up along with the post-shock settling (see, e.g., 
Mayle et al. 1987, Myra & Burrows 1990, Bruenn 1993, 
Bruenn et al. 1995). 

Just like the neutrino emission from the less opaque, 
hot mantle region of the protoneutron star dominates the 
supernova neutrino fluxes for several ten to some hundred 
milliseconds after core bounce, more than 90% of the neu- 
trino emission of the merger comes from the extended, 
hot torus and only a minor fraction of less than 10% orig- 
inates from the very opaque and dense central core of the 
merged object. Figures ||, ||, |[ and || display the local 
emission rates of v e , v ei the sum of u x , and the sum of 
all neutrino types in the orbital plane for models A64, 
A128, B64, and C64, respectively, at a time when quasi- 
stationary conditions have been established. Figures and 

give the corresponding information in two orthogonal 
cut planes vertical to the equatorial plane for the repre- 
sentative models A128 and C64. One can clearly see that 
the surface-near regions, in particular the toroidal cloud 
of matter surrounding the central, compact body, emits 
neutrinos at much higher rates. 

Using the values from these figures, one can esti- 
mate the relative importance of core and disk emission. 
The dense core regions (p > 10 13 g/cm 3 ) inside a ra- 
dius of R c ~ 20 km lose neutrino energy with a typi- 
cal rate of Q £ ~ 10 32 erg/cm 3 /s and account for a lu- 
minosity of about L„ ;C w (47r/3)it! 3 (2e w 3.4 • 10 51 erg/s, 
while the emission from the surrounding torus-shaped 
disk region with outer radius i?d « 40 km, height 2h rj 
40 km, and typical total neutrino energy loss rate Q £ ~ 
5 • 10 32 erg/cm 3 /s emits a luminosity of roughly L u < \ w 
{2nRlh - 47ri? 3 /3)Q £ w 8.4 • 10 52 erg/s. The sum of 
and I/„ : d is of the size of the results displayed in Fig. [l| 

Figures |5|— |l~0| also demonstrate that the disk emits v e 
fluxes that are larger than the v e fluxes. Electron neu- 
trinos are most strongly emitted from surface-near re- 
gions where the optical depth to v e by absorption on neu- 
trons [y e + n — > p + e~) in the neutron-rich matter is 
smallest. In contrast, the production of electron antineu- 
trinos by positron captures on neutrons and of heavy- 
lepton neutrinos via electron-positron pair annihilation 
(e _ + e + — > v + D) requires the presence of large numbers 
of positrons and therefore occurs predominantly in those 
parts of the disk which have been heated most strongly 
and have a low electron degeneracy. Since v e are absorbed 



on the less abundant protons in the inverse /3 + process and 
the opacity to v x is primarily caused by neutrino-nucleon 
scatterings only, electron antineutrinos and heavy-lepton 
neutrinos can escape on average from deeper and hotter 
layers than electron neutrinos (see Figs. p[-^0|). This ex- 
plains the higher mean energies of the emitted v e and v x . 

As expected from the very similar neutrino luminosi- 
ties and nearly equal mean energies of emitted neutrinos, 
the neutrino emission maps do not reveal major differences 
between the models A64, B64, and C64. Most neutrino 
emission comes from the disk region where the neutrino 
optical depths are lower, and in Figs. ||, ||, and ^| one can 
clearly recognize the high-density, very opaque inner part 
of the core of the merger from its roughly two orders of 
magnitude smaller energy loss rates. The cuts perpendic- 
ular to the orbital plane (see Fig. |Io|) show the ring-like 
main emission region that surrounds the central core and 
has a banana- or dumb-bell shaped cross section. Model 
A128 has significantly more fine structure but the overall 
features and characteristics of the neutrino emission do 
not change with the much better numerical resolution of 
this simulation. 

3.2. Thermodynamics and composition 

Figure [ll] displays contour levels in the temperature- 
density plane of the electron neutrino chemical potential 
(measured in MeV) for an electron fraction of Y e = 0.046. 
In Fig. [l^ the contours corresponding to vanishing elec- 
tron neutrino chemical potential /i„ e = are plotted for 
different values of Y e . Figures [ll] and |l2| provide informa- 
tion about how chemical equilibrium is shifted according 
to the equation of state of Lattimer & Swesty (1991) when 
neutron star matter changes its temperature-density state. 

Cold neutron star matter at neutrinoless /3-equilibrium 
is in a state with fi Ue — 0. When such gas of the out- 
ermost ~ 0.1-0.2 Mq of the neutron star with density 
p 10 14 g/cm 3 and electron fraction Y e <, 0.025 (see 
Fig. 2 in Paper I) is expanded and heated while the lep- 
ton fraction Y\ cp = Y e + Y Ve « Y e stays roughly constant 
as it is the case for a fast change where neutrino losses 
are too slow to compete, the /3-equilibrium is shifted into 
the region of negative /i„ e values. This usually implies 
that the electron degeneracy of the matter is drastically 
decreased, too, because the electron chemical potential 
H e = jU We + ji n — fi p also drops when n Ve attains negative 
values. For hot gas (i.e., gas at conditions lying above the 
nose- like feature of the curves in Figs, [ll] and |l^) the state 
with /i„ e = corresponds to a higher value of Y e which is 
seen in Fig. [l2] by moving along lines parallel to the or- 
dinate. On the lepton-number loss time scale associated 
with neutrino emission, the gas will tend to evolve again 
towards the /3-equilibrium with \x Vs = by an enhanced 
production and emission of D e relative to v e . Notice that 
from Figs. |l] and one infers that the same arguments 
are true for the case that neutron star matter with den- 
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Fig. 11. Contours of the neutrino chemical potential (mea- 
sured in MeV) in the temperature-density plane for fixed elec- 
tron fraction Y e . The contours are spaced with increments of 
10 and are labeled with their respective values 



Fig. 12. Contours of vanishing neutrino chemical potential 
= 0) in the temperature-density plane for different values 
of the electron fraction Y e as indicated by the line labels 



sities above p « 10 14 g/cm 3 is strongly compressed. From 
the properties of the high-density equation of state we 
therefore deduce that during the coalescence of neutron 
star binaries the hot gas in the compact, compressed core 
region of the merger as well as the heated, decompressed 
disk matter will radiate u e more copiously than v e . 

This explains the relative sizes of electron neutrino and 
antineutrino luminosities as discussed in Sect. 3.1. Driven 



by this imbalance of the emission of v e and v e , the initially 
very neutron-rich matter (Y e ^ 0.095 everywhere in the 
neutron star, see Fig. 2 in Paper I) gains electron lepton 
number and becomes more proton-rich again. Figure p"3| 
shows this evolution for model A128 from the start of the 
simulation until its end at 8.80 ms. The snapshots of the 
Y e distribution in the orbital plane visualize how, as a con- 
sequence of the rapid neutrino loss from the disk region, 
Y e in this region climbs from initial values of 0.02-0.06 to 
values of more than 0.18 in some parts. In the core re- 
gion the neutrino emission proceeds much more slowly so 
that Y e changes only slightly during the simulated time. If 
we continued our computations for a long enough time to 
see the matter cooling again by neutrino losses (provided 
the configuration is stable for a sufficiently long period), 
this process of Y e increase would again be inverted and 
the gas would evolve towards the cold, deleptonized, very 
neutron-rich state again. 

Figure |lj displays the final situations (t ss 10-11 ms) 
in the models A64, B64, and C64 to be compared with 
panel f of Fig. O. As in case of the neutrino emission, one 



notices very similar properties of all four models. In model 
A64 the peak Y e values in the disk region are as high as 
0.22. Note that in all four models very neutron-rich mat- 
ter is swept off the grid, a tiny fractio n o f which might 
potentially become unbound (see Sect. 2.3 and Paper I). 
It is also interesting to see the still very elongated and de- 
formed neutron-rich inner region of models A64 and B64 
which indicates ongoing strong dynamical and pulsational 
activity of the massive core. This is not so pronounced in 
model C64 where the anti-spin setup of the initial model 
has caused the dissipation of a large fraction of the rota- 
tional energy during the coalescence of the neutron stars. 
Model A128 has also a much more circular core region be- 
cause, as described in Paper I, the better resolution allows 
for a much more fine-granular flow pattern which contains 
a large fraction of the initial vorticity and kinetic energy 
in small vortex structures. 

Figure |l5| presents a collection of plots of parameters 
that give information about the thermodynamical state 
in the orbital plane of model A128 and about the nu- 
clear composition of the gas at the end of the simulation. 
Panels a and b show the electron degeneracy parameter 
Ve — Pe/(kBT) (fcs is the Boltzmann constant, /i e the 
electron chemical potential) and the electron neutrino de- 
generacy parameter r\ Vii = ^/(fceT), respectively. Con- 
cordant with the discussion above, the /3-equilibrium con- 
ditions in the whole star are characterized by r\ Vs 
(panel b). In the core values between —3 and —6 can be 
found, while in the disk moderately negative values are 
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Fie. 13. Time evolution of the SDatial distribution of the electron fraction K, in the orbital Diane of model A128. The times 
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Fig. 14. The spatial distributions of the electron fraction 
Y e in the orbital plane for models A64, B64, and C64 at the 
end of the calculations (times given in the top right corners). 
The contours are linearly spaced with intervalls of 0.02 and the 
bold contours are labeled with their respective values. The grey 
shading emphasizes the levels with darker shading indicating 
higher values of Y e . The plots have to be compared with panel 
f of Fig. [[3] which shows the same information for model A128 
at time t = 8.80 ms 



present (around —1) and in some regions the medium has 
evolved back to a state close to r\ Va w 0. Comparison with 
panels a and b of Fig. || shows that in these regions the 
emission rates of v e and v e are already very similar again 
whereas the production of electron antineutrinos is clearly 
dominant in those parts of the disk with the most nega- 
tive values of r] Ve . The electron degeneracy is moderate 
(r\ e ps 2-3) in the disk but climbs to numbers around 
r\ e w 25 near the center. 

The temperature k^T and entropy are displayed in 
panels c and d, respectively, of Fig. [lf| Detailed infor- 
mation about the evolution of the temperature in all 
models was given in Paper I (Figs. 4-7, 14-17, and 20, 
21). At the end of the simulation model A128 has the 
highest temperatures of k-gT ps 30 MeV in hot spots 



located in a shell around the central high-density core 
where k-gT w 10 MeV. The disk has been heated up 
to k-gT w 10 MeV in regions of density p sa 10 12 - 
10 13 g/cm 3 , k B T > 6 MeV where p £ 10 n g/cm 3 , and 
/cbT <; 1-2 MeV for p *J 10 10 g/cm 3 . The corresponding 
entropies are less than 1 fce/nucleon in the core region 
and between 3 and slightly more than 7 fcs/ nucleon in 
the disk (panel d) . In model C64 similar disk entropies are 
found while in model A64 specific entropies up to about 
9 fce/nucleon and in model B64 up to even 10 fce/nucleon 
develop towards the end of the simulated evolution. 

Such high entropies allow only minor contributions of 
nuclei to be present in the gas in nuclear statistical equi- 
librium at the densities found for the disk matter on our 
computational grid. Most of the nuclei are completely dis- 
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Fig. 15. Cuts in the orbital plane of model A128 showing different thermodynamical and composition parameters at the end 
of the simulated evolution (time t = 8.80 ms). Panel a gives the electron degeneracy parameter r; e (contours linearly spaced 
with steps of one unit), panel b shows the electron neutrino degeneracy parameter rj Vc (contours linearly spaced with steps of 
0.5 units), panel c displays the temperature distribution (contours linearly spaced with steps of 2 MeV), panel d is a plot of 
the entropy per nucleon (contours linearly spaced with steps of 1 Ub, /nucleon) , and panel e informs about the mass fraction of 
a particles in the medium (contours logarithmically spaced in steps of one dex) 
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integrated into free nucleons (the mass fraction of heavy 
nuclei is below the lower limit of ~ 10 -8 returned from 
the equation of state of Lattimer & Swesty (1991)), and 
only small admixtures of a particles are possible. Panel e 
of Fig. 15 shows that the mass fraction X a of a particles 
is typically less than about 10 -3 . Only in the outermost 
parts of the disk X a ~ 10 -2 because there the tempera- 
tures are low enough, k-gT ~ 1-2 MeV, that some of the 
free nucleons can recombinc. 



4. Neutrino-antineutrino annihilation 

Neutrino-antineutrino annihilation in the surroundings of 
the merger has been proposed to create a sufficiently ener- 
getic fireball of e + e _ -pairs and photons to explain gamma- 
ray bursts at cosmological distances. We attempt to put 
this idea to a quantitative test. With the given informa- 
tion about the fluxes and spectra of the neutrino emission 
of all grid cells (see Paper I for technical details), it is 
possible to evaluate our hydrodynamical models for the 
energy deposition by ^^-annihilation in a post-processing 
step. Since the neutrino luminosities become large only 
after the merging of the two neutron stars and in partic- 
ular after the gas torus around the compact central body 
has formed, we consider the late stages of our simulated 
merger evolutions as the most interesting ones to perform 
the analyses. In the phase when quasi-stationary condi- 
tions have been established, the neutrino luminosities have 
reached their saturation levels and the annihilation rates 
have become maximal. 



4--1- Numerical evaluation 

Neglecting phase space blocking effects in the phase spaces 
of e~ and e + , the local energy deposition rate (energy 
cm _3 s _1 ) at a position r by annihilation of Vi and z?j 
into e + e _ -pairs (which is the dominant reaction between 
neutrinos and antineutrinos) can be written in terms of 
the neutrino and antineutrino phase space distribution 
functions f V( = f Vi (e, n, r, t) and f Si = f 0i (e',n',v,t) 
as (Goodman et al. 1987, Cooperstein et al. 1987, Janka 
1991) 
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When the energy integrations are absorbed into (energy- 
integrated) neutrino intensities I Vi and I Pi , 
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The integrals over fl and f2' sum up neutrino and 
antineutrino radiation incident from all directions. 9 
is the angle between neutrino and antineutrino beams 
and (e) Ui and (e) Pi are suitably defined average spec- 
tral energies of neutrinos and antineutrinos, respectively. 
The weak interaction cross section is a = 1.76 x 
10 -44 cm 2 , m e c 2 = 0.511 MeV is the electron rest- 
mass energy, c the speed of light, and the weak cou- 
pling constants are (C% + Ci) v p 



(c v + c A y 



2.34, C : 



= (Cv ~ CaT + 
\(2Cl-C\) « 1.06 



and (d + d) Vm9m = (C v - CaT + (C v + C A - 2) 2 » 
0.50, C Wx = | [2(C V - l) 2 - (C A ~ I) 2 ] » -0.16 for 
v x = Vpi v T and C A — \ , Cy = i + 2 sin 2 6>w with 
sin 2 #vv = 0.23. The total energy deposition rate at the 
position r is given as the sum of the contributions from 
annihilation of v e and v e , and 9^, and v T and v T : 



(4) 



When working with a discrete grid the integrals in 
Eq. (||) are replaced by sums over all cells k, 



b dn I v — > ^2 A ^fc ' I ",t 

in i 



(5) 



AJ7fe is the solid angle with which cell k is seen from a 
position r at distance <4 = |r — | when r*. is the loca- 
tion of the center of cell k. In order to avoid the need to 
take into account projection effects, we define an effective 
radius D associated with the cells of the cartesian grid 
used for the hydrodynamical modelling by setting the cell 
volume 14 = AxAyAz = (Ax) 3 equal to the volume of a 
sphere V = 4ttL> 3 /3: 



1/3 

D = [ f - I Ax 

Air , 



With the projected area A = 7rD 2 we obtain 



Afi fc 



-in 



2/3 



Ax 



(6) 



(7) 
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neutrino species vi = u e) u e , v^, v^, 
I Vi ,k is therefore given by 



v T or v T , the intensity 



1/3 



Eq. (0) and Eq. (§) allow us to evaluate the sum of Eq. 
sufficiently accurately by 



4 



(9) 



With the average energy of neutrinos emitted from cell 
k, ( e )vi,k, (see Paper I) and the angle enclosed by the ra- 
diation from cells k and k' at the position r, cos dkw = 
(r - v k ) (r - r k >) / (d k d k >), the integrals of Eq. (§) finally 
become sums over all combinations of grid cells k with 
all cells k' . This double sum has to be evaluated at 
all positions r where the energy deposition rate by vv- 
annihilation is to be determined. The corresponding com- 
putational load is appreciable but can be significantly re- 
duced by taking into account only those grid cells which 
emit towards r, i.e. whose emission into the outward di- 
rected half space around n p is also pointing to position r. 
The criterion for this being fulfilled is n p • (r — r^) > 0. 

4-2. Numerical results 

The numerical post-processing procedure described in the 
previous section yields the total energy deposition rate 
Qvv by neutrino- antincutrino annihilation as a function 
of the position r. Figure |l^ shows the results, averaged 
over azimuthal angles, in a quadrant of the r-z-plane per- 
pendicular to the orbital plane for the three models A64, 
B64, and C64. The evaluation was performed only in that 
region around the central part of the merger where the 
local energy loss rate by neutrino emission is less than 
10 30 erg/cm 3 /s and the density is below 10 11 g/cm 3 . Den- 
sity levels are indicated by dashed contour lines in the 
plots. 

One can see that in all three models the highest rate 
of energy deposition (> 10 30 erg/cm 3 /s) occurs in the 
outer regions of the disk (35 km < r < 45 km in the 
orbital plane) within about 25 km above and below the 
orbital plane. Here the energy deposition rate is much 
larger than the energy loss rate. Because this deposition 
transfers energy into gas layers with densities of still more 
than 10 10 g/cm 3 , one must suspect that a baryonic wind 
will be created very similar to the neutrino-driven wind 



co 

CD 



Using the simplifying assumption that a grid cell radi- 
ates neutrinos with isotropic intensity into the half space 
around the outward direction defined by the local density 
gradient n p = Vp/|V/o|, the flux j Vi ,k is related to the 
neutrino radiation intensity I Vi ,k by j Vi> k = 7r ^ i ,fc- With 
an effective emissivity Q^(i>i) (see Paper I) which repre- 
sents the energy emission of cell k per cm _3 s _1 in a single 
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Fig. 17. Cumulative total energy deposition rates by neu- 
trino-antineutrino annihilation into electron-positron pairs for 
the three models A64, B64, and C64 as functions of the 
equatorial distance r from the center of the merger once 
for the case that the integration is done for polar distances 
20km < |z| < 64km and another time for 0km < \z\ < 64km. 
The integrals were performed over the values of the local en- 
ergy deposition rates plotted in the three panels of Fig. Il6] 



caused by neutrino energy deposition in the surface lay- 
ers of the nascent neutron star in a type-II supernova (for 
information about the neutrino-driven wind from form- 
ing neutron stars, see Duncan et al. 1986, Woosley & 
Baron 1992, Witti et al. 1994, Woosley et al. 1994; see 
also Hernanz et al. 1994). Instead of creating a fireball of 
a nearly pure relativistic electron-positron-photon plasma 
which might lead to an energetic gamma-ray burst, this 
energy is used up to lift baryons in the strong gravita- 
tional field of the massive central body. Consequently, the 
expansion velocities of this matter are nonrelativistic. If 
too much of this wind material gets mixed into the pair- 
photon plasma, the baryonic load M w can become too 
high to allow for Lorentz factors T in the required range 
E u9 /{M w c 2 ) £ 100. 



of T = 1/V1 ~{v/c) 2 

The large energy deposition rates at radii between 
30 km and 50 km in the equatorial plane and at moderate 
heights above and below the orbital plane can be explained 
by the closeness to the main neutrino emitting ring in the 
disk between 20 km and 30 km (see Figs. 1-0). Accord- 
ing to Eqs. (||), (||), and (g) the annihilation rate decreases 
at least with the fourth power of the distance to the neu- 
trino radiating grid cells k and kl: Q„ 9 oc l/{d\d\,). Be- 
cause of the influence of the geometrical factor (1 — cos 0) 2 
in Eq. (j|) (this factor accounts for the dependence of the 
annihilation probability on the relative velocity of the in- 
teracting neutrinos and the kinematically allowed phase 
space for the reaction), Q^ p decreases even more steeply 
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when one moves away from the source and large-6* colli- 
sions between neutrinos become less and less frequent. At 
distances from the merged object d = \Jr 2 + z 2 <; 50- 
60 km we find that along each radial beam the local en- 
ergy deposition rate as a function of d becomes a power 
law with power law index between —6 and —8. Towards 
the rotation (z-) axis, the growing distance to the neu- 
trino producing disk region cannot be compensated by 
the higher chance of head-on collisions between neutrinos. 
This is the reason why the contour levels in the plots of 
Fig. |l^ bend towards the polar region of the merger near 
the z-axis. 

The integral values of the energy deposition rate by 
^^-annihilation in the surroundings of the merger out to 
equatorial distances r are shown in Fig. [IT]. For each of the 



annihilation rates, B64 
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Fig. 16. Maps of the total local energy deposition rates 
(in erg/cm 3 /s) by neutrino-antineutrino annihilation into elec- 
tron-positron pairs in the vicinity of the merger for the three 
models A64, B64, and C64 at time t = 10 ms after the start 
of the simulations. In one quadrant of the r-z plane orthogo- 
nal to the orbital plane (at z = 0), the plots show values ob- 
tained as averages of the energy deposition rate over azimuthal 
angles. The corresponding solid contour lines are logarithmi- 
cally spaced in steps of 0.5 dex, the grey shading emphasizes 
the levels with dark grey meaning high energy deposition rate. 
The dashed contours indicate levels of the azimuthally aver- 
aged density, also logarithmically spaced with intervals of 0.5 
dex. The energy deposition rate was evaluated only in that 
region around the merged object, where the mass density is 
below 10 11 g/cm 3 and the energy loss rate by neutrino emis- 
sion is smaller than 10 30 erg/cm 3 /s. One can see that due to 
the closeness to the main neutrino radiating disk region, most 
annihilation energy is deposited in the "outer" (in r-direction) 
and "upper" (in ^-direction) regions of the disk (dark grey ar- 
eas) 

models A64, B64, and C64 we have evaluated the spatial 
integral once for vertical heights 20 km < \z\ < 64 km and 
another time for km < \z\ < 64 km. The maximum up- 
per integration limit is determined by the largest distances 
where the local rates Q^p were calculated. However, one 
can see from Fig. [l7] that the curves tend to approach a 
saturation level for r — > 64 km from which we conclude 
that extending the integrations into the region d > 64 km 
would not change the results significantly. From a com- 
parison of both cases one can recognize that only a minor 
fraction of about 1/5-1/4 of the annihilation energy is 
converted into e + e~-pairs in the region above and below 
the disk. Only at heights \z\ <; 20 km is the baryon density 
low enough that most of the converted energy might end 
up in a relativistic fireball. However, in the whole region 
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km < \z\ < 64 km the energy deposition rate is rather 
small, only about 2-4 • 10 50 erg/s; the "useful" fraction is 
probably less than 1/4 of that. The models were evaluated 
at times when the neutrino emission of the models had al- 
ready achieved a maximum value and a quasi-stationary 
state. At earlier times the neutrino luminosities are much 
lower and therefore the integral values of the energy depo- 
sition rate are even smaller than those displayed in Fig. [l7[ 
For the disk or torus geometry of our models the an- 
nihilation efficiency, defined as e V y = {AE vj> /6.t)/{L Vc + 
Lv e + 4L„ M ) where AE V[ ,/At is the total energy deposition 
rate by w-annihilation, can be found to be of the order 
of a few tenths of a per cent: 



(2. ..3) • 10" 



20 km 



1.5 • 10 52 erg/s 13MeV R d 



(10) 



The factors in this equation can be deduced by simple an- 
alytical transformations and dimensional considerations of 
the volume integral of the annihilation rate (Eq. (|l0|)). R^ 
is the inner radius of the disk, the radius which roughly 
separates the compact inner core from the more dilute 
gaseous cloud of matter around. The inserted numbers are 
representative values taken from Figs. |l|, ||, and [|. The ef- 
ficiency of Eq. (^) for converting emitted neutrino energy 
into electron-positron pairs by w-annihilation is in good 
agreement with what is obtained in supernova simulations 
or analytical estimates assuming emission in a spherical 
geometry instead of the disk-dominated emission of our 
numerical models. 

From the results discussed here we conclude that there 
is no chance to obtain the energies needed for a cosmologi- 
cal gamma-ray burst by ^-annihilation during the instant 
of the merging of binary neutron stars. Before the hot 
cloud of gas around the central, dense body has formed, 
the neutrino luminosities stay rather low and only very 
little energy is deposited by ^-annihilation. Even later 
when a disk around the compact, merged body has formed 
and the neutrino emission has reached a high level, the 
energy deposition rates of 2-4 • 10 50 erg/s lead to a total 
deposited energy of only a few 10 48 erg within the com- 
puted time of 10 ms. One would need the strong neutrino 
emission to be maintained for periods of about 10 s to 
pump an energy of more than 10 51 erg into a fireball of 
e + e~-pairs and photons. Note that these statements are 
not changed if the observed gamma-ray bursts are beamed 
events. If one relies on the emission geometry of our mod- 
els, a constrained solid angle Af2 into which the fireball ex- 
pands and becomes visible in a jet-like outburst would also 
(and by the same factor) reduce the spatial volume where 
^-annihilation in the surroundings of the merger could 
deposit the useful energy. Correspondingly, the numbers 
given in Fig. [Tt] would have to be scaled down by a factor 
Af2/ (4?r). In view of the fact discussed in Sect. 2.3 that the 
central, compact core of the merger has a mass that can 
neither be stabilized by internal pressure nor by rotation 
and should therefore collapse into a black hole on a time 



scale of milliseconds, we conclude that the energy available 
by pair annihilation of neutrinos is lower than the desired 
and required "canonical" value of 10 51 /(47r) erg/steradian 
by more than a factor of 1000. 

So far we have considered the phase of the merging of 
binary neutron stars and the evolution that follows im- 
mediately afterwards. This, however, leaves the question 
unanswered whether a disk around the most likely forming 
black hole could produce and emit neutrinos on a much 
longer time scale and could thus power a gamma-ray burst 
by the discussed mechanism of neutrino-pair annihilation? 
Since simulations of this scenario are beyond the present 
capabilities of the employed numerical code and, in addi- 
tion, the results depend on the unknown viscosity in the 
disk, we attempt to develop a simple model of the behavior 
and properties of such a disk with respect to its neutrino 
emission and the strength of ^-annihilation. 



4-3. Simple model for the post-merging emission from the 
disk 

Our simulations suggest that some material, possibly 
about 0.1 M , could remain in a disk around the central 
black hole. This disk will be heated by viscous dissipation 
and will emit neutrinos and antineutrinos until its matter 
is accreted into the central black hole. The efficiency of 
w-anmhilation increases linearily with the v luminosity 
(Eq. JlO|)) and thus a short, very luminous neutrino burst 
is more efficient to create an e + e~-pair fireball than the 
same energy emitted on a longer time scale with smaller 
neutrino fluxes. It has to be investigated whether enough 
energy can be provided in the pair-photon fireball by the 
neutrino emission from the disk to explain a 7-ray burst at 
cosmological distances. Viscosity effects have a crucial in- 
fluence on the disk evolution and on the neutrino emission. 
Viscous forces, on the one hand, transfer angular momen- 
tum between adjacent fluid elements and determine the 
accretion time scale and accretion rate. Viscous dissipa- 
tion of rotational energy, on the other hand, heats the 
disk and is thus essential for the neutrino emission. Disk 
size, disk temperature, disk viscosity, and neutrino emis- 
sion properties can therefore not be chosen independent of 
the accretor mass and disk mass. In the following we shall 
attempt to relate these quantities by simple considerations 
and conservation arguments. 

The lifetime of the disk will decrease with larger dy- 
namic viscosity 77 because the viscous force that gener- 
ates a torque carrying angular momentum outward is in- 
creased. For a (Newtonian) Keplerian disk the viscous 
force (per unit area) in the angular (0) direction, is 
simply expressed by the component t r § of the viscous 
stress tensor (see, e.g., Shapiro & Teukolsky 1983): 



-t r 4> — - 7?f2k , 



(11) 
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where fik = ^/GM/r 3, is the Keplerian angular velocity. 
The torque T exerted by the viscous stress is given by 
T = fyr (2-kt 2h) when 2h ~ 2R S is taken as the ver- 
tical diameter of the thick disk around the black hole 
with mass M and Schwarzschild radius R s — 2GM/c 2 w 
9(M/3M Q )km. The accretion rate dM/dt can be esti- 
mated (roughly) by setting the viscous torque equal to 
the rate dJ/dt ~ (dM / dt)Okr 2 at which angular momen- 
tum is consumed by the black hole due to the accretion of 
matter from the disk: 



dM 
~df 



6irr/R s 



(12) 



From that, the accretion time scale of a disk of mass AAfd 
is estimated to be 



AMd 
6iri]R s 



(13) 



Thus, the lifetime of the disk is determined by the out- 
ward transport of angular momentum through the vis- 
cous torque. Equation ( |l3| ) shows that it decreases with 
the value of the dynamic viscosity as I/77. 

Viscous dissipation generates heat in the disk at a rate 
per unit volume of (see Shapiro & Teukolsky 1983) 



dQ 
dt 



(14) 



At steady-state conditions the maximum dissipation rate 



occurs at a radius r m 1.36 i?d 



• Rd when Rd is the 



inner radius of the disk which is taken to be the innermost 
stable circular orbit around the central accreting black 
hole, R d w 3R S w 6GM/c 2 . Using this in Eq. @) one 
obtains for the maximum rate at which frictional heat is 
liberated, 



dQ 
dt 



1 GM 
28 V ~rT 



(15) 



Thus, the viscous heating rate increases linearly with r\. 

For small viscosity 77 the viscous heating time scale 
is long and the disk remains rather cool, also because 
cool matter is comparatively transparent for neutrinos and 
therefore the neutrino cooling time scale is short. In that 
case the neutrino luminosity for a disk with volume Vd is 
L v ~ Vd dQ/dt oc 77 and the total energy radiated in neu- 
trinos, E v ~ L v t acc , becomes independent of 77 because of 
iacc oc I/77. With Eq. (|ic| ) one finds that the energy con- 
verted into e + e~ by ^-annihilation increases proportional 
to 77: E V y = e vi) E v cx L^iacc oc 77. In the optically thin case 
the mean energy of emitted neutrinos, (e„), which enters 
the calculation of e„ p will also increase with 77 and cause 
a slightly steeper than linear dependence of E vj> on 77. 

If 77 is large, the disk is heated rapidly and strongly 
and thus becomes opaque for neutrinos. With a neu- 
trino diffusion time scale idiff that is much longer than 
the heating time scale the neutrino luminosity is L v ~ 



Vd (dQ/d£)/j a ccAdiff which is only indirectly dependent 
on 77 through idiff and thus the (viscosity dependent) 
gas temperature T. In that case E v oc l/(rftdis) and 
E u0 oc l/(?7^ iff ). Note that the average energy of emitted 
neutrinos, (e„), which also determines e^p, is only very 
weakly dependent on the viscosity of the disk in the op- 
tically thick case because it reflects the conditions at the 
neutrino decoupling sphere (see Eq. ( |2l| ) below). The dif- 
fusion time scale increases with the disk temperature and 
thus with the disk viscosity due to the energy dependence 
of the weak interaction cross sections. This leads to a de- 
crease of E v! > with 77 that is steeper than I/77. 

The considerations above suggest that the annihilation 
energy E vj> has a pronounced maximum at a particular 
value 77* of the dynamic viscosity. Because of E v0 cx L 2 
the annihilation of neutrinos and antineutrinos is more ef- 
ficient when a certain energy is emitted in a short time 
with a high luminosity rather than over a long period 
with a moderate flux. If 77 is small, L v stays low. If 77 
is very large and the interior of the disk very hot and thus 
neutrino-opaque, the neutrino luminosity L v scales with 
the inverse of the neutrino diffusion time scale and with 
the total energy E v that can be emitted in neutrinos dur- 
ing the lifetime t acc of the disk. This energy E v decreases 
in case of very large 77 because t acc becomes shorter and 
the internal energy cannot be completely radiated away in 
neutrinos before the neutrino-opaque matter is accreted 
into the black hole. The kinetic energy that is converted 
into internal energy by viscous friction is entirely trans- 
ported away by neutrinos and the fluxes are largest, if the 
diffusion time scale is similar to the accretion time scale 
but not much longer. The optimum value 77* is therefore 
determined by the condition t acc idiff- 

Let us assume that the part of the disk where most of 
the neutrinos are emitted has a mass AMd and is a ho- 
mogeneous torus with center at 4i? s and radius R s (inner 
radius 3i? s , outer radius 5R S ) (Mochkovitch et al. 1993; 
Jaroszyriski 1993). This is a fairly good picture in view of 
the shape and structure of the disk that we obtained in our 
numerical simulations?]. In terms of R s the volume of the 
disk torus is Vd = S-pR^ and its surface Sd = 167r 2 i? 2 . 
With the neutrino mean free path A = p<T c s/m u and 
AMd = pVd the diffusion time scale is approximately 
given by 



idiff 



cA 



3(AM d )a cff 
8ir 2 R s m u c 



(16) 



Here m u is the atomic mass unit and the thermally av- 
eraged effective neutrino interaction cross section <j c h = 
Y) . YiGi is defined as the sum of the cross sections <7j times 



1 In the simple model considered here, the structure and 
geometry of the disk is assumed to be given. It is not self- 
consistently determined in dependence of the gas temperature 
and thus in dependence of the competing effects of viscous 
heating and neutrino cooling. 



24 



M. Ruffert et al.: Coalescing neutron stars - a step towards physical models 



the number fractions Yi — rii /ub of the corresponding re- 
action targets for all neutrino processes in the medium, 
i.e., neutrino scattering off n, p, e~ , e + , charged-current 
absorptions of v e and v e by n and p, respectively, and 
impair interactions. We find 

2 



Ccff 



(1...4) • 10" 



k B T 
5MeV 



cm 



(17) 



the exact value depending on the neutrino type, the neu- 
trino degeneracy and neutrino spectra, and the detailed 
composition of the medium. For the entropies, densities, 
and temperatures obtained in our simulations the gas in 
the disk is completely disintegrated into free nucleons; nu- 
clcon as well as lepton degeneracy plays a negligible role 
(see Sect. 3.2). Therefore fermion phase space blocking ef- 



fects are unimportant. The thermal average of the neu- 
trino cross section was evaluated by using a Fermi-Dirac 
distribution function with a vanishing neutrino chemical 
potential, p v = 0. In case of incomplete dissociation of 
the nuclei the neutrino opacity should still be within the 
uncertainty range associated with the cross section varia- 
tion of Eq. (H). Setting t acc (Eq. @) with AM d = pV d 
equal to tdiff (Eq. Jl6|)), one determines the value of the 
shear viscosity in the disk, where ^-annihilation yields 
the largest energy, as 

-2 



4-7T m u c 



9 (JcB 



(1.7...6.9)-10 2 



k B T 
5MeV 



.(18) 



The range of values accounts for the uncertainty in the 
effective neutrino interaction cross section. For the typical 
composition of the disk material the cross section is more 
likely near the upper limit of the given interval, in which 
case the lower value of rf is favored. 

Let us now consider a disk with this optimum value rf . 
Making use of rf , the interior temperature of the disk can 
be estimated by setting the integral rate of viscous energy 
generation in the disk, L visc ~ Vd(dQ/dt), equal to the 
luminosity due to neutrino diffusion, L v ~ V d E v /t d {tf ~ 
(^c\)S d (e L ,/R s ). Here £„/i? s is an approximation to the 
gradient of the neutrino energy density in the disk and 
£ v ~ 3 • | tiradTJnt * s * ne sum °f * ne ener gy densities of all 
three kinds of non-degenerate w-pairs with a ra( j being the 
radiation constant. One finds 



5.6(AM d , i) 1/4 -R i 



3/4 
s,9 



5MeV 

where AM d)0 i = AM d /(0.1M o ) and i? s , 9 = i? s /(9km) 
is normalized to the Schwarzschild radius of a 3 Mq black 
hole. Note that due to the dependence of £ d iff and of rj* on 
(Toff (according to Eqs. (16) and ([L8|), respectively) k B T- mt 
does not depend on the neutrino interaction cross section 
and is therefore insensitive to its uncertainty. 

Plugging the result of Eq. ( |l9| ) for the interior disk 
temperature into Eq. (|l^) yields for the optimum disk vis- 
cosity 



V 



(5.5...22)-10 25 (AAf d ,oi)- 1/2 i? s 3/ 9 2 



cms 



(20) 



the interval of values again corresponding to the range of 
possible values of the effective neutrino interaction cross 
section, rf from Eq. (|2(]) can now be used in Eq. (15) 
to calculate L v isc ~ V d (dQ/di), which, when set equal 
to the neutrino luminosity expressed in terms of temper- 
ature and surface area S d of the neutrinosphere, L v ~ 
S , d (3|a rad T s 4 urf )c/4, leads to an estimate of the neutri- 
nospheric temperature 



5McV 



(0.7...1.0)-(AM d ,oi)- 1/8 i? s 1/ 9 8 



(21) 



The temperature of the neutrino emitting disk surface is 
around 5 MeV and rather insensitive to the exact value 
of the effective neutrino interaction cross section (slightly 
larger result for smaller cross section), to the disk mass 
M d , and to the inner disk radius R d ~ 3R S . 

The optimum value rf for the dynamic viscosity as 
given in Eq. ( pp| ) corresponds to an effective a-parameter 
of a w r]*/(pc s R s ) - (1.7... 7.0) • 10" 3 when R s = 
9 km, AM d = 0.1 Mq, and p = AM d /{8n 2 R*) = 3.44 • 
10 12 g/cm 3 are used, and the sound speed c s is evaluated 
with k B T = 28 MeV. For these values of density and tem- 
perature the gas pressure is dominated by relativistic par- 
ticles, i.e., photons, electrons, positrons, and neutrinos. 
Neutrino shear viscosity does not contribute significantly 
to rf . In the neutrino-opaque case it is estimated to be 



1 e„\ _ e u m u 

3 C 3/9(T c ffC 



(1 ...4) • 10 



2:i 



k B T y 
5 MeV/ 



10 12 g/cm 3 \ g 



P 



cm s 



(22) 



where the interval of the numerical value is again asso- 
ciated with the uncertainty of the effective cross section 
CT ofT . For p = 3.44 • 10 12 g/cm 3 and k B T = 28 MeV one 
finds rj y ~ (0.9... 3.6) • 10 24 gcm _1 s _1 . A temperature as 
high as k B T ~ 70-80 MeV is required for the neutrino 
viscosity to become large enough to account for rf . 

For the diffusion and accretion time scales one obtains 
by inserting Eq. (pC|) into Eq. 



t dm - (53 ... 212) (AM d , i) 3/2 i? s . 9 5/2 ms . (23) 



This time is much longer than the dynamical time scale 
(0(1 ms)) and the neutrino equilibration time scale (^ 
lms). Therefore our assumptions that neutrinos diffuse 
in the disk and are in equilibrium with the matter are 
confirmed a posteriori. The total neutrino luminosity is 



(0.62 ... 2.48) • 10 53 (AA/ d ,oir 1/2 i? s 5/ 9 2 erg/s . (24) 



In Eq. ( |23| ) the smaller values and in Eq. ( pi] ) the larger 
ones correspond to the case of larger viscosity rf and thus 
smaller neutrino cross section cr c ff according to the pos- 
tulated equality of Eqs. ( |l3| ) and (|l6|). The total energy 
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E v = L v t acc radiated away over the time t acc is indepen- 
dent of both and becomes 



7 3i? s - 42 (AMd)c 



1.3 ■ 10 52 AM d , i erg 



(25) 



which is (approximately) equal to the Newtonian gravita- 
tional binding energy EWd = ^GM(AM d )/R d = \E glav + 
E ro t | of mass AMj at the inner disk radius Rd = 3i? s 
where the matter is swallowed by the black hole (E grav 
is the gravitational potential energy E mt the rotational 
energy). Here it is assumed that no rotational kinetic en- 
ergy is extracted from the black hole which is equivalent 
to a zero stress boundary condition at Rd- This requires 
that within R d the gas spirals into the black hole rapidly 
without radiating, an idealization which is probably justi- 
fied (see, e.g., Shapiro & Teukolsky 1983). The small dis- 
crepancy between the factors i and j of our calculation 
results from the fact that we consider a simple one-zone 
model of a homogeneous disk. We find that the radia- 
tion efficiency of the disk in our simplified treatment is 
£^/(AM d c 2 ) = |j « 7.5% (exact value for a thin, New- 
tonian accretion disk: ss 8.3%). This result has to be 
compared with the radiation efficiency of about 5.7% for 
relativistic disk accretion onto a nonrotating black hole 
and with the radiation efficiency of 42.3% for a maximally 
rotating black hole with a prograde accretion disk (see 
Shapiro & Teukolsky 1983). Since our numerical simula- 
tions suggest the formation of a central black hole with 
a relativistic rotation parameter a = Jc/{GM) that is 
clearly less than 1 (Ruffert et al. 1996), the reference value 
for the radiation efficiency for disk accretion onto a non- 
rotating black hole is relevant and our Newtonian disk 
evolution model most likely overestimates the amount of 
energy that can be carried away by neutrinos before the 
accreted mass finally plunges rapidly from Rd to the event 
horizon. 

Using the results of Eqs. (|2l]) and (|24| ) to compute the 
^-annihilation efficiency e„p according to Eq. (1^) and 
employing the integral energy E v emitted in neutrinos as 
given from Eq. ( |25| ) , one can obtain a result for the energy 
Euv = e V vE v deposited in an e + e _ -pair-photon fireball by 
the annihilation of v and v radiated from the disk. With 
R d = 3i? s = 27 km, (e Ve ) w 3fc B r sU rf ~ (10.5... 15.0) MeV, 
and L Ve ~ (| ... i) L u (see Sect. 3T), i.e, L Ve (1 . 
10 52 erg/s, we get 



•4) 



E VB ~ (1.1...9.4)-10 49 (AM d , 01 ) 3/8 i? s 13 9 /8 erg 



(26) 



The upper and lower bounds of the interval for E V y corre- 
spond to the most extreme (maximum and minimum, re- 
spectively) choices for e„ p , L Ve , and (e„ e ) . Notice that the 
analytical estimates of the neutrino luminosity (Eq. (pi[)) 
and the mean energy of neutrinos emitted from the disk, 
(e„) fa 3fc B T surf w 11. ..15 MeV (Eq. (pTh ) , agree well with 



our numerical results for the phase shortly after the merg- 
ing. While the dynamical time scale of the merging is of 
the order of 1 ms and the post-merging evolution was fol- 
lowed by our numerical simulations for a period of about 
10 ms, the disk emits neutrinos with similar luminosities 
for a much longer time of a few hundred milliseconds 
(cf. Eq. (p3|)). Therefore, Eq. ( p6| ) gives a number for the 
energy deposition by ^-annihilation that is a factor of 
10-100 larger than the ~ 10 48 erg calculated in Sect. 4.2 



5. Discussion 

In this paper we have reported about hydrodynamical 
calculations of the merging of equal-mass binary neu- 
tron stars with different initial spins. We have anal- 
ysed the models for their neutrino emission, for neutrino- 
antincutrino annihilation in the surroundings of the 
merger, and for the thermodynamical conditions in the 
merged object. 

5.1. Mass loss and nucleosynthesis 

The dynamical merging proceeds within a few milliseconds 
after the simulations were started from an initial center- 
to-center distance of 42 km of the two 1.6 Mq neutron 
stars. Shortly after the two stars have fused into one com- 
pact object with a mass of about 3 Mq and an average 
density of more than 10 14 g/cm 3 , spun-off matter forms a 
less dense toroidal cloud (p w 10 12 g/cm ) that is heated 
to temperatures of 5-10 MeV by friction in shock waves 
and strong pressure waves sent into the surrounding gas by 
the oscillations and periodic pulsations of the central high- 
density body. Roughly 0.1 Mq of material receive a large 
momentum and are pushed beyond the grid boundaries. 
However, only a small fraction of at most 10 _5 ...10 -3 Mq 
of the matter has a total energy (internal plus kinetic plus 
gravitational) large enough to allow the gas to become 
unbound. In model A128 which has the best numerical 
resolution, in fact none of the matter can escape from the 
gravitational potential of the merger, in contrast to model 
A64 where it is several 10 -5 Mq. Notice, however, that en- 
ergy released by the recombination of free nucleons into 
nuclei — an effect which is taken into account in our sim- 
ulations by the use of the "realistic" equation of state of 
Lattimer & Swesty (1991) — and energy production by 
nuclear reactions proceeding in the cooling and expand- 
ing gas (Davies et al. 1994) can aid the mass ejection and 
could increase the unbound mass relative to our estimates. 

The ejection of matter is also very sensitive to the 
amount of angular momentum that is present in the merg- 
ing binary. The largest mass ejection was found in model 
B64 immediately after the merging, because the initial 
configuration of this model had a solid-body type rota- 
tion and thus the largest specific angular momentum in 
regions far away from the system axis. In contrast, model 
C64 showed the largest mass loss a few milliseconds later. 
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The initial anti-spin setup of this model led to vigorous 
vibrations of the central body and to the outward accelera- 
tion of material some time after the two neutron stars had 
formed a single object. Whether and how much mass can 
be dynamically lost during the post-merging evolution, 
however, will also depend on the stability of the merged 
object. The central, compact core of the merger is so mas- 
sive that it can be stabilized neither by internal pressure 
for the currently favored supranuclear equations of state, 
nor by its rapid rotation (for details, see Sect. 4.1.3 of 
Paper I) . One therefore has to expect its collapse into a 
black hole within a few milliseconds. In case of model C64 
not much matter would be expelled if the gravitational in- 
stability sets in before the large-amplitude post-merging 
oscillations have taken place. 

The less dense cloud of gas that surrounds the mas- 
sive, very dense central body is stabilized by internal pres- 
sure because its rotational velocities are significantly less 
than the Kepler velocity (see Fig. 12 of Paper I). When 
the massive core collapses into a black hole, the ambient 
matter will therefore be swallowed up by the black hole on 
a dynamical time scale of 10 -4 ...1CP 3 s. Only gas with a 
sufficiently large angular momentum will have a chance to 
remain in a disk or extended torus around the black hole. 
From there it will spiral into the black hole on the much 
longer time scale of viscous angular momentum transport. 
One can estimate that with typical orbital velocities of 
about 0.25c as found in our models, only gas at radii be- 
yond about 44-107 km has an angular momentum that is 
large enough (for more information, see Sect. 4.1.3 of Pa- 
per I) . This minimum orbital radius is outside of the grid 
boundaries of our simulations. We therefore conclude that 
essentially all the mass on the computational grid will dis- 
appear in a forming black hole more or less immediately, 
and only the ~ 0.015-0.15 Mq of material that have been 
swept off the grid might be able to end up in a toroidal 
"disk" around the central black hole. A disk mass of about 
0.1 Mq should be taken as an extreme upper limit for the 
considered scenario. Like the dynamically ejected mass, 
the amount of gas that ends up in a disk is sensitive to 
the initial angular momentum of the neutron star binary 
and to the relative times of black hole formation and mass 
spin-off. Last but not least, a quantitative answer seems 
to depend also on the numerical resolution, with the trend 
that better resolved models yield smaller estimates for the 
possible disk mass and mass ejection. 

A mass of 10~ 4 ...10~ 3 Mq that is dynamically ejected 
during the merging of binary neutron stars might have 
important implications for nucleosynthesis (Lattimer & 
Schramm 1974, 1976; Eichler et al. 1989). Dependent on 
the phase when the mass loss occurs, the expelled gas will 
start its expansion from different initial conditions of en- 
tropy and composition. The ejection of initially very cool, 
low-entropy material might be caused by the tidal inter- 
action during the last stages of the inspiral and during the 
mass transfer phase of very close non-equal mass binaries. 
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Fig. 18. Expected average mass numbers of nuclei formed by 
an r-process starting with NSE conditions for low entropies and 
with the a-process for high entropies. The plot shows contours 
of constant average mass number in the Y e -s plane (s in ks 
per nucleon). Y e and s define the thermodynamical conditions 
in the expanding gas at a temperature of about 5 • 10 9 K, i.e., 
before the r-processing takes place. The dynamical evolution is 
characterized by the time scale r of the adiabatic expansion be- 
tween T ~ 7- 10 9 K and T ~ 3- 10 9 K. Solid lines correspond to 
a time scale of 50 ms, dashed lines to 100 ms. Between the con- 
tours for mass numbers A = 120 and A = 210 the shaded area 
marks the region where a suitably mass-weighted combination 
of the r-process yields for different conditions will produce a 
solar-system like abundance pattern. For larger values of Y e , 
no significant r-processing can take place, for lower values of 
Y e very strong r-processing will primarily lead to nuclei in the 
region of the actinides 



If the two components have "nearly equal" initial masses, 
the mass transfer is unstable and within a few orbital pe- 
riods the lighter star can be completely dissipated into a 
thick, axially symmetric disk around the primary. Some 
fraction of the surface material might escape the system 
(Lattimer & Schramm 1974, 1976). If the initial mass ra- 
tio of the two stars is large, the binary is stable against 
dynamical-time scale mass transfer and there is the in- 
teresting possibility that the secondary (the less massive 
component) is stripped to the minimum mass of stable 
neutron stars, at which stage it will explode (Page 1982; 
Blinnikov et al. 1984, 1990; Eichler et al. 1989; Colpi et 
al. 1989, 1991, 1993; Colpi & Rasio 1994). However, recent 
investigations suggest that stable mass transfer is unlikely 
because the initial mass of the secondary must already 
be very small (below ~ 0.4 M Q ; Bildsten & Cutler 1992, 
Kochanek 1992, Rasio & Shapiro 1994, Lai et al. 1994). In 
addition, an unreasonably high value of the neutron star 
viscosity is needed to enforce corotation and to maintain 
tidal locking, because stable mass transfer requires a dy- 
namically stable Roche limit configuration which can only 
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exist in synchronized systems with extreme mass ratios. 
Such systems are essentially ruled out for neutron stars 
(Lai et al. 1994). If mass shedding in the discussed situa- 
tions occurs, it would lead to the ejection of initially cold, 
very low-entropy and very neutron-rich material. Subse- 
quent radioactive /3-decays of unstable, neutron-rich heavy 
nuclei that are present in the decompressed matter will 
heat the expanding gas to temperatures around 0.1 MeV, 
which will give rise to r-process conditions (Lattimer et 
al. 1977, Meyer 1989, Eichler et al. 1989). 

Gas spun off the exterior parts of the dilute toroidal 
cloud that surrounds the compact core of the merged 
binary, instead, has been heated to temperatures of *S 
1 MeV by friction during the merging and post-merging 
evolution. This heating has produced entropies of a few 
fcs per nucleon. Neutrino emission has already raised Y e 
from initial values Y e ss 0.02 to slightly less neutron-rich 
conditions with Y e ss 0.05-0.2. The electron degeneracy 
is only moderate, rj e « 2. All these parameters are very 
similar to the conditions found in the shocked outer lay- 
ers of the collapsed stellar core in a type-II supernova 
where the site of the classical r-process has been suggested 
(see, e.g., Hillebrandt 1978). Compared with the super- 
nova case, the range of F e -values in the potentially ejected 
merger material is on the low side. An r-processing oc- 
curring under such neutron-rich conditions would be very 
efficient and should preferentially produce r-process nu- 
clei with very high mass numbers (Fig. |l^). The predomi- 
nant production of high-mass r-process elements would be 
in concordance with the fact that one cannot expect the 
formation of all Galactic r-process material in the con- 
sidered low-entropy ejecta. With our numerical estimates 
of a few 10 -4 Mq for the mass loss per merger event and 
with the possible event rate of NS-NS and NS-BH mergers 
of 10 -6 ...10 -4 per year per galaxy (Narayan et al. 1991, 
Phinney 1991, Tutukov et al. 1992 Tutukov & Yungelson 
1993), which corresponds to about 10 4 ...10 6 events during 
the lifetime of the Galaxy, only 1-100 Mq of the ~ 10 5 M Q 
of Galactic r-process elements could be produced. But 
if the r-process were strong enough, all the Galactic ac- 
tinides (e.g., about 40-50 Mq of Th) might be accounted 
for by the material shed during neutron star merging! 

However, the neutron-rich wind that is driven by neu- 
trino energy deposition in the outer disk regions can also 
contribute to the nucleosynthetic input into the interstel- 
lar medium. With a total neutrino luminosity of about 
10 53 erg/s the wind will have a mass outflow rate of the or- 
der of 0.001-0.01 Mqs" 1 (Qian & Woosley 1996, Woosley 
1993b), depending on the gravitating mass of the merger, 
the disk mass and geometry, and the neutrino emission as 
a function of time. For a duration of the outflow between 
some fractions of a second and a few seconds, one might 
therefore have another 10~ 4 Mq up to several 10 -2 Mq of 
material that are expelled with very interesting thermo- 
dynamical properties. 



Like the neutrino-driven wind from new-born neutron 
stars, the neutrino-heated material should have signifi- 
cantly higher entropies than the matter that is dynami- 
cally ejected from the disk or torus by momentum transfer 
during core pulsations. Owing to the particularities of ge- 
ometry, gravitational potential, and neutrino emission in 
the merger situation, the expansion time scales as well as 
the degree of neutronization might be significantly differ- 
ent from the supernova case. Since the neutrino emission 
from the core and the disk of the merger is dominated 
by the v e fluxes, absorptions of 9 e in the wind material 
(p e + p — > n + e + ) will be more frequent than the ab- 
sorption of less abundant v e {v e + n — > p + e~) and will 
keep the expanding material neutron-rich. Because of a 
larger luminosity ratio L Pc /L l , i! « 3-4.5 but similar mean 
v e and v e energies with (ep e ) / (e Ve ) ~ 1.5-1.8, the expand- 
ing wind will have a lower electron fraction than in the 
supernova case (see Qian & Woosley 1996 for a discus- 
sion of the electron fraction in neutrino-driven winds). 
Values as low as Y e w [1 + Lp e (e Pe )/(_L„ e (e„ e ))] « 0.1- 
0.2 seem possible in the neutrino wind from the merger. 
For such low values of Y e a strong r-processing can hap- 
pen even at modest entropies of s ~ 50-100 /ce/nucleon 



(see Sect. 5.2 for an estimate of the wind entropies) which 
are too low to allow for the formation of r-process nu- 
clei in the neutrino-driven winds from protoneutron stars 
for the typical electron fractions of Y e 0.35-0.4 found 
there (Witti et al. 1994, Takahashi et al. 1994, Woosley et 
al. 1994, Qian & Woosley 1996). Figure [l8] visualizes this 
and shows that even for rather slow expansions with ex- 
pansion time scales of more than 100 ms nuclei with mass 
numbers A between 150 and 210 can be formed. 

Unfortunately our current models allow only rough es- 
timates of the mass loss and the conditions in the ejected 
matter. Our simulations could not directly follow the ejec- 
tion of mass from the merger because they suffered from 
the limitations due to the use of the computational grid 
and due to an insufficient numerical resolution, especially 
of matter at low densities. A detailed and meaningful 
analysis of the very interesting aspects of a possible r- 
processing in the dynamically ejected low-entropy mate- 
rial and in the high-entropy neutrino-driven winds from 
NS-NS or NS-BH mergers has to be postponed until mod- 
els are available which yield more quantitative informa- 
tion about the long-time evolution of the merged object, 
the torus geometry, and the neutrino-matter interactions 
in the outer parts of the torus. Only such models can 
give evidence about the duration of the mass loss and the 
amount of material that is ejected with different entropies, 
different expansion time scales, and different degrees of 
neutronization, all of which determine the nucleosynthetic 
processes (see Fig. |l| and also Woosley & Hoffman 1992, 
Witti et al. 1994). 
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5.2. Neutrino emission and gamma-ray bursts 

The luminosities and mean energies of the neutrinos emit- 
ted from merging neutron stars are very similar to those 
calculated for supernovae and protoneutron stars. After 
the two neutron stars have merged, luminosities up to sev- 
eral 10 52 erg/s are reached for every neutrino species and 
the average energies of v e leaking out of the merger are 10- 
13 MeV, of v & they are 19-21 MeV, and of heavy-lepton 
neutrinos around 26-28 MeV. However, the neutrino emis- 
sion exhibits characteristic differences from the supernova 
case, too. 

The total neutrino luminosity from merging neutron 
stars does not increase to a value above 10 52 erg/s before 
the hot, toroidal gas cloud around the dense and com- 
pact core of the merger starts to form. More than 90% 
of the peak neutrino emission of about 10 53 erg/s stems 
from the "disk" region where the optical depths and thus 
the neutrino diffusion time scales are significantly smaller 
than in the core. Another pecularity is the fact that the 
very neutron-rich, decompressed and heated neutron star 
matter predominantly emits electron antineutrinos. This 
should hold on until the electron fraction in the medium 
has grown to a level where the increase of the lcpton num- 
ber by v e emission is compensated by the v e losses. If 
the merged configuration remained stable for a sufficiently 
long time, the deneutronization phase will be superseded 
by an extended period where the heated gas deleptonizes 
and cools again and thus evolves back to the state of 
cold neutron star matter. However, it is very likely that 
the merged object, which contains essentially the baryonic 
mass of two typical neutron stars (about 3Mq), does not 
remain gravitationally stable. For all currently favored nu- 
clear equations of state it should collapse to a black hole 
long before the cooling is finished. 

If the gravitational instability of the massive core of 
the merged object sets in before the hot gaseous torus has 
formed and if all the surrounding gas falls into the black 
hole immediately, the neutrino emission from the merger 
will stay fairly low with a total luminosity of less than 
r-j 10 52 erg/s. This is much too low to get sufficient energy 
for a cosmological gamma-ray burst from ^^-annihilation 
during the final stages of the inspiral of the two neutron 
stars and during the first 1-4 ms right after the merging. 
The annihilation efficiency of neutrinos and antineutrinos, 
e V y — (dEvp/dt)/ L„, increases proportional to the neu- 
trino luminosity (Eq. (|Io|)) and the neutrino energy de- 
position rate E ue with the product of neutrino and an- 
tineutrino luminosities. At the time of maximum neutrino 
emission some 6-8 ms after the stars have merged, we 
calculate an annihilation efficiency of (2-3) • 10 -3 and a 
neutrino energy deposition rate of 2-4 • 10 50 erg/s in the 
whole space outside the high-density regions of the com- 
pact core and the surrounding "disk" . The integrated en- 
ergy deposition during the simulated evolution of about 
10 ms is therefore less than 4 • 10 48 erg (assuming maxi- 



mum neutrino fluxes during the whole considered times). 
Even with the unrealistic assumption that all the vv- 
annihilation energy could be useful to power a relativis- 
tic pair-photon fireball, this energy would fail to account 
for the canonical ~ 10 51 /(47r) erg/steradian of a typical 
gamma-ray burst at cosmological distances (e.g., Woods 
& Loeb 1994; Quashnock 1996) by nearly three orders of 
magnitude. Of course, these arguments are not conclusive 
if there is strong focussing of the expanding fireball into 
a narrow solid angle <5f2. In this case an observer would 
deduce a largely overestimated value (by a factor Att/SVI) 
for the energy in the fireball and in the gamma-ray burst 
if he assumed isotropy of the emission. However, for the 
considered merger scenario and the geometry of the post- 
merging configurations in our simulations, it is very hard 
to imagine how the required strong beaming of the fireball 
into a jet-like outflow could be achieved. 

It is interesting to note that if the neutrino emission 
calculated for our merger models would continue for a few 
seconds, which is a typical duration of observed gamma- 
ray bursts (e.g., Norris et al. 1994, Kouveliotou 1995), a 
burst energy of about 10 51 erg could well be accounted for 
by the annihilation of neutrinos and antineutrinos. An ac- 
cretion disk or torus around the central black hole could 
provide a luminous neutrino source for the required period 
of time. This time span is much longer than the times cov- 
ered by our hydrodynamical modelling. Since the current 
numerical simulations were neither able to mimic the ef- 
fects of a central black hole nor to track the evolution of 
the merger for a sufficiently long time, we attempted to 
develop a simple analytical model in Sect. |4.3| to give us 
insight into the principal dependences of the energy depo- 
sition by the annihilation of neutrinos emitted from the 
disk or torus. 

This torus model was based on Newtonian physics and 
did not determine the torus structure self-consistently. 
However, the analytic treatment took into account the 
effects of viscous angular momentum transport, viscous 
heating, neutrino cooling, and partial neutrino opaque- 
ness. To first order and on a qualitative level, our consider- 
ations should also be valid for accretion disks around black 
holes in general relativity (see, e.g., Shapiro & Teukol- 
sky 1983) and for tori around Schwarzschild black holes 
(see, e.g., Chakrabarti 1996). In fact, the employed as- 
sumptions about the torus geometry are supported by 
general relativistic investigations of neutron tori (Witt 
et al. 1994; Jaroszynski 1993, 1996) and the neutrino lu- 
minosities from our simple torus model are compatible 
with those obtained from the relativistic analyses. More- 
over, the neutrino flux and neutrinospheric temperature 
calculated analytically are also in good agreement with 
the results of our numerical models for the post-merging 
phase when the neutrino emission has reached its satura- 
tion level. 

Equation (|2^) gives the estimate of the ^^-annihilation 
energy E v „ for our analytical disk model. We find that E vi > 
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could at best lie between 1.1 • 10 49 erg and 9.4 • 10 49 erg 
for a disk of 0.1 M Q around a 3M© black hole. The life- 
time of such a disk is determined by the time scale of the 
outward transport of angular momentum and is estimated 
to be several ten up to a few hundred milliseconds. The 
ranges of values account for the uncertainties in the neu- 
trino opacity (which depends on the composition of the 
medium and on the neutrino spectra) and for the corre- 
sponding variation of the neutrino luminosity and mean 
energy of the emitted neutrinos. Further uncertainties due 
to the unknown disk viscosity are circumvented by deriv- 
ing the result for a value ij* of the dynamic shear viscosity 
(Eq. (^0j)) which assures a maximum result for E v „. 

Unless focussing or beaming of the expanding pair- 
photon fireball towards the observer plays an important 
role, the annihilation energy of Eq. ( p6| ) is too low by 
more than a factor of 10 to explain powerful cosmolog- 
ical gamma-ray bursts. Weak bursts with an energy of 
about 10 50 erg and durations of less than or around a sec- 
ond, however, do not seem to be completely excluded on 
grounds of Eqs. (|6|) and (|23|). The result of Eq. (|2^) and 
in particular the upper value of ~ 10 50 erg should be con- 
sidered as a very optimistic maximum estimate for the 
energy deposited by ^-annihilation in the surroundings 
of the disk. Equation ( p6| ) was derived for the most favor- 
able conditions and by making a whole sequence of most 
extreme assumptions, the combination of all of which ap- 
pears rather unlikely 

In the first place it was assumed that the dynamic vis- 
cosity of the disk adopts the optimum value rj* of Eq. fll8|). 
For much smaller viscosities the toroidal disk should stay 
rather cool and the neutrino fluxes correspondingly low. 
For much larger values of the viscosity the lifetime of the 
disk will decrease because the mass accretion rate into the 
black hole increases with the rate of the outward transport 
of angular momentum mediated by viscous forces. In this 
case the neutrino luminosities will be bounded by the fact 
that the viscous friction will heat up the disk to very high 
temperatures and thus the neutrino absorption and scat- 
tering cross sections, which scale roughly with the square 
of the neutrino energy, will increase. Therefore the neu- 
trino diffusion time scale will increase, too, and the neu- 
trino cooling will become inefficient. As a consequence, 
most of the dissipated gravitational and rotational energy 
could be advected into the black hole when the matter, 
after having lost part of its angular momentum, spirals 
in through the innermost stable circular orbit (advection- 
dominated regime). 

Equation (^6|) also represents an optimistically high 
value of the energy deposition by ^-annihilation because 
the radiation efficiency of about 8% obtained for our New- 
tonian model of a Keplerian accretion disk is an upper 
bound to the radiation efficiency of a relativistic disk 
around a nonrotating black hole where it is less than 6%. 
Moreover, it should be remembered that the annihilation 
efficiency e^p of Eq. (|l0|) most likely overestimates the 



useful fraction of the annihilation energy by a consider- 
able factor. Our hydrodynamical simulations show that 
the fraction of the neutrino energy deposited in the pos- 
sibly baryon-poor region above and below the disk but 
not in the plane of the disk (where it will serve to drive 
a baryonic, nonrelativistic wind instead of creating a rel- 
ativistically expanding pair-photon fireball) could be as 
small as some 20-25% of the number given in Eq. (26) 
(cf. Fig. [It]). Even more, general relativistic effects were 
neglected in the numerical simulations and analytical con- 
siderations presented here. Jaroszyhski (1993) showed that 
they lower the energy that can be transported to infinity 
significantly (by about 80%). 

For the whole uncertainty range of the neutrino inter- 
action cross section (Eq. fll7|)) and for the corresponding 
range of neutrino luminosities and mean energies of emit- 
ted neutrinos, we find that the result of Eq. (^6|) falls short 
of the desired value -E„ p /(47r) ~ 10 51 /(47r) erg/steradian 
by at least an order of magnitude if the disk mass is of the 
order of 0.1 M and the central black hole has a mass of 
about 3 M Q . From our numerical models of NS-NS merg- 
ers, one concludes that a disk with a mass close to or even 
larger than 0.1 M might be formed only under very spe- 
cial conditions. Even if there is a high angular momentum 
in the system due to neutron star spins as in our model 
B64, the amount of material that has a chance to form a 
disk is hardly as much as 0.1 M (see Paper I). In models 
A64 and C64 the lower angular momentum allows only 
little material to possibly remain in a disk. Only matter 
outside of the boundary of our computational grid has 
enough angular momentum to be rotationally stabilized. 
Before some matter in models A64 and C64 has acquired 
a sufficiently large angular momentum to be lost off the 
computational grid — a process that is partly aided by 
pressure waves created by the wobbling and ringing of the 
central, compact object — the merger, however, has prob- 
ably already collapsed to a black hole which swallows up 
most of the surrounding, pressure supported matter on a 
dynamical time scale. 

Massive, self-gravitating tori around black holes may 
be subject to general relativistic global instabilities that 
lead to catastrophic runaway mass loss and may provide 
by far the shortest evolutionary time scale of such tori as 
recently argued for stationary polytropic tori by Nishida 
et al. (1996) and for stationary neutron tori by Nishida 
& Eriguchi (1996). Such an instability would have the 
same implications as the case of extremely large disk vis- 
cosity discussed above where the accretion time scale of 
the torus due to the rapid viscous angular momentum 
transport could be much smaller than the neutrino dif- 
fusion time scale. As a consequence of the rapid accretion 
of the torus material, most of the internal energy of the 
gas would be carried into the black hole along with the 
gas instead of being radiated away by neutrinos. Because 
the duration of the neutrino emission would be very short 
without a compensating increase of the neutrino luminos- 
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ity, the total energy emitted in neutrinos would be much 
smaller than in the extreme and optimum situation con- 
sidered in the derivation of Eq. (E6I) in Sect. 4.3. There- 



fore, global disk instabilities might be another threat to 
neutrino-powered gamma-ray bursts. However, it has still 
to be demonstrated whether global runaway instabilities 
develop in the non-stationary situation and how they be- 
have and evolve in the presence of changes of the angular 
momentum distribution. 

Even if the lifetime of the merger and of the neu- 
trino radiating accretion torus is long enough and neu- 
trino annihilation could provide a powerful "engine" for 
creating a fireball, yet another major concern for the vi- 
ability of the considered 7-ray burster scenario comes 
from the baryonic wind that is blown off the surface of 
the merger and accretion torus by neutrino heating. This 
neutrino-driven wind is unavoidable when large neutrino 
fluxes are emitted and a small fraction of these neutri- 
nos annihilate or react with nucleons in the low-density 
gas in the outer layers of the merger. In order to ob- 
tain bulk Lorentz factors T ~ ^/(MyC 2 ) 100 one 
can only allow for mass loss rates M w ^ 2 • 10~ 6 Af s _1 
if the rate at which the pair-photon fireball is supplied 
with energy is E uD « 3 • 10 50 erg/s. With about 2/3 or 
~ 2 • 10 50 erg/s of this energy being transferred to the 
dilute out er regions of the accretion torus in our models 
(Sect. 4.2 ) by j/jz-annihilation (neglecting additional heat- 
ing by neutrino-electron scattering and neutrino absorp- 
tion), we compute a mass loss rate of at least M w ^ ^E v p ■ 
[GAf/(3i? s )] _1 £ 7-10- 4 A/ Q s- 1 for a black hole with mass 
M = 3 M Q and Schwarzschild radius R s ss 27 km. From 
this we get E vD / '(M w c 2 ) Ss 0.25 and estimate an entropy 
s w (e + P) I '{ks,T p 1 'm u ) (s internal energy density, P pres- 
sure) of s sa [4m u c 2 /(3fc B T)] [(£W/3)/(M w c 2 )] £ 100- 

200 fc B /nucleon when fc B T ~ 0.5 ... 1 MeV. Therefore, un- 
less the neutrino-driven wind can be hindered to penetrate 
into the pair-photon fireball, e.g., in a region along the sys- 
tem axis by centrifugal forces, there is no chance to obtain 
highly relativistic fireballs with T <; 100. 

These issues might also be critical when a neutron star 
merges with a black hole or when two non-equal mass neu- 
tron stars coalesce. Simulations indicate (Lee & Kluzniak 
1995) that already after the dynamical interaction of a 
neutron star and a black hole a cloud of baryonic material 
might "pollute" the surroundings even near the system 
axis. Moreover, instead of an accretion torus, a stable bi- 
nary system might form with a more massive black hole 
circulated by a low-mass neutron star companio n w hich 
is possibly unstable to explosion (see also Sect. 5A and 
Blinnikov et al. 1984, Eichler et al. 1989). Future, better 
resolved computations covering a wider range of equations 
of state and masses of the interacting stars will have to 
show the conclusiveness of their results. 

It is interesting to speculate whether there is a chance 
to get tori more massive than in the merging of two 



(nearly) equal-mass neutron stars when a neutron star 
merges with a black hole or when a small neutron star 
is tidally disrupted before merging with a companion 
neutron star which has a significantly larger mass (e.g., 
Eichler et al. 1989, Narayan et al. 1992, Mochkovitch et 
al. 1993, and references therein). On the one hand, a more 
massive disk can radiate neutrinos for a longer time than 
the accretion time scale of the innermost, most strongly 
neutrino radiating part near the last stable orbit around 
the black hole. With much more matter being at larger 
radii, the neutrino emitting torus region considered and 
normalized to a mass of 0.1 Mq in the derivation of 
Eq. (^6|) will be continuously refed by material advected 
inward. On the other hand, Eq. (^6|) suggests that the en- 
ergy that can be provided in the pair-photon fireball by 
^^-annihilation increases steeply with the mass M of the 
central black hole, E vP oc Rl 3 ^ 8 oc Af 13 / 8 . Larger gravitat- 
ing masses at the center, e.g., a massive stellar black hole, 
might therefore allow for more powerful 7-ray bursts, at 
least, if one relies on the simplified picture developed in 
the derivation of Eq. ( p6| ) which does not take into account 
the dynamics of the tidal interaction between neutron star 
and black hole and its implications for the disk formation 
and the disk mass. 

For the reasons outlined above it appears to be ex- 
tremely hard to account for the energies of cosmological 
7-ray bursts by ^-annihilation, at least in the case of 
merging binary neutron stars if beaming or focussing of 
the fireball towards the observer does not play a crucial 
role. The same conclusion was arrived at by Jaroszyhski 
(1993, 1996) who investigated models of relativistic tori 
around rotating stellar mass black holes and tested dif- 
ferent values of specific angular momentum, viscosity, and 
entropy. He found that the neutrino emission and annihila- 
tion energy from these tori is insufficient to explain the en- 
ergies of cosmological 7-ray bursts except for tori around 
Kerr black holes with very high angular momentum, i.e., 
with relativistic rotation parameters a ~ 1. This might in- 
deed be realized in the collapsed cores of rapidly rotating 
Wolf-Rayet stars in the "failed supernova" or "collapsar" 
scenario (Woosley 1993a). However, models for the stel- 
lar evolution of Wolf-Rayet stars seem to indicate that 
the angular momentum loss during the mass-loss phases 
is too large to allow for the formation of very rapidly 
spinning black holes (Langer, personal communication). 
Moreover, with ^-annihilation as the source of the en- 
ergy of the pair-photon fireball, these models might come 
into additional trouble if, as claimed by Quashnock (1996), 
the homogeneity of the distribution of the bursts in the 
BATSE 3B Catalog and the non-association of the bursts 
with large-scale structures of luminous matter in our ex- 
tragalactic neighbourhood implies such large distances of 
the 7-ray burst sources that the required total energy out- 
put in 7-rays is larger than 3 • 10 51 erg. 

Our investigations do not include effects due to pos- 
sible convective overturn and instabilities in the disk 
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that might occur as a result of specific entropy (or com- 
position) inversions caused by neutrino effects or vis- 
cous heating. Such dynamical processes were recently 
found to be present in multi-dimensional models of type- 
II supernovae (Herant et al. 1994; Burrows et al. 1995; 
Janka & Muller 1995, 1996) and were computed for 
axisymmetric advection-dominated accretion flows with 
two-dimensional hydrodynamical simulations by Igumen- 
shchev et al. (1996) (sec also Chen 1996). Although the 
disk remains globally stable, shorter-wavelength modes 
may affect the flow dynamics and effective disk viscos- 
ity significantly. In addition, such instabilities might in- 
crease the neutrino fluxes and the average neutrino en- 
ergies considerably and thus might help ^^-annihilation. 
Our three-dimensional hydrodynamical simulations were 
not followed for a sufficiently long time to see whether 
such overturn processes occur in the merging scenario, 
and the simplified one-zone torus model does not take 
into account an enhancement of the neutrino fluxes by 
possible convcctivc transport. Also, magnetic fields in the 
merging neutron stars and in the torus were disregarded. 
Within a lifetime of a few tenths of a second, initial in- 
fields of ~ 10 12 -10 13 G might be amplified by a factor of 
100 or more in the rapidly rotating disk around the black 
hole (rotation periods ~ 1 ms) and might become energet- 
ically important (Rees, personal communication). A rel- 
ativistic magneto-hydrodynamical wind of extremly high 
luminosity, perhaps associated with a binary neutron star 
merger, was suggested to generate 7-ray bursts by Thomp- 
son (1994). In this respect NS-BH mergers and neutron 
star collisions and coalescence need further theoretical 
investigation. Moreover, neutron star collisions were re- 
cently pointed out as potential origin of short cosmological 
7-ray bursts by Katz & Canel (1995a) who argued that hot 
matter fragments could be ejected and, when they become 
neutrino-transparent on dynamical time scales, could lead 
to the emission of neutrinos with very high luminosities. 

6. Summary and conclusions 

The neutrino emission and neutrino-antineutrino anni- 
hilation during the coalescence of binary neutron stars 
were investigated. To this end the three-dimensional New- 
tonian equations of hydrodynamics were integrated by 
the Riemann-solver based "Piecewise Parabolic Method" 
on an equidistant Cartesian grid with a resolution of 
64 x 64 x 32 or 128 x 128 x 64 zones. The properties of neu- 
tron star matter were described by the equation of state 
of Lattimer & Swesty (1991). Energy loss and changes of 
the electron abundance due to the emission of neutrinos 
were taken into account by an elaborate "neutrino leakage 
scheme" . We have simulated the coalescence of two iden- 
tical, cool (initially k^,T c w 7MeV) neutron stars with a 
baryonic mass of about 1.6 M©, a radius of 15 km, and an 
initial ccnter-to-center distance of 42 km for three differ- 
ent cases of initial neutron star spins. 



The total neutrino luminosity prior to and during 
the dynamical phase of the coalescence is very small 
(L u & 10 51 erg/s), becomes about 1-2 • 10 52 erg/s when 
the stars have merged into one rapidly spinning massive 
body, and climbs to 1-1.5 • 10 53 erg/s after spun off ma- 
terial has formed a hot toroidal cloud with a mass of 
0.1-0.2 Mq around the wobbling and pulsating central 
object. The neutrino fluxes are clearly dominated (~90- 
95%) by the emission from this "disk". Since the disk 
matter is neutron-rich, v e are radiated with a luminosity 
that is a factor 3-6 higher than the (individual) luminosi- 
ties of v e and v x (= v^, v^, v T , v T ). The mean energies 
of the emitted neutrinos are very similar to those of su- 
pernova neutrinos, (e„ e ) w 12MeV, (e Se ) « 20MeV, and 
27 MeV. 

When the neutrino luminosities arc highest, only about 
0.2-0.3% of the energy emitted in neutrinos is deposited 
in the immediate neighborhood of the merger by vv- 
annihilation, and the maximum integral energy deposition 
rate is found to be about 3-4 • 10 50 erg/s. Thus, to pump 
an energy of the order of 10 51 /(47r) erg/steradian into a 
fireball of e + e~-pairs and photons, the strong neutrino 
emission would have to continue for several seconds. Since 
a collapse of the central core of the merger with a mass 
of <; 3 Mq into a black hole within milliseconds seems 
unavoidable, we conclude that the available energy is in- 
sufficient by a factor of about 1000 to explain gamma-ray 
bursts at cosmological distances. However, it appears pos- 
sible that an accretion torus with a mass of <~ 0.1-0.2 M© 
remains around the central black hole and is accreted on 
the time scale of viscous angular momentum transport. 
Analytical estimates suggest that even under the most 
favorable conditions in this torus and with an optimum 
value of the disk viscosity, annihilation of vv pairs emitted 
from this torus provides an energy that is still more than 
a factor of 10 too small to account for powerful cosmo- 
logical gamma-ray bursts, unless focussing of the fireball 
expansion plays an important role. 

A few 10~ 4 Mq of very neutron-rich, low-entropy mat- 
ter may be dynamically ejected shortly after the neutron 
stars have merged, and another 10~ 4 up to a few 10~ 2 M© 
of strongly neutronized, high-entropy material might be 
carried away from the accretion torus in a neutrino-driven 
wind on a time scale between a fraction of a second and 
a few seconds. The contamination with these baryons is 
a severe threat to a relativistic fireball. Aspects of nu- 
cleosynthesis in these ejecta were discussed. Because of 
the neutron-richness of the ejected material and the dom- 
inance of the v e luminosity from the merged object and 
its accretion torus, conditions suitable for the formation 
of r-process elements might be realized more easily than 
in the neutrino wind from newly formed neutron stars. 

It seems to be very difficult to fulfill the energetic re- 
quirements of cosmological gamma-ray bursts with the an- 
nihilation of vv pairs emitted from an accretion disk or 
torus around a stellar mass black hole. If ^-annihilation 
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is nevertheless to be saved as energy source for relativistic 
pair-photon fireballs — despite of the problems exposed 
by our numerical and analytical results and the critical 
issues addressed in the discussion of Sect. |— then one is 
forced to consider the following possibilities. 

The neutrino luminosities from the accretion torus 
could be considerably higher than obtained in our mod- 
els, but the mechanism to achieve this has yet to be iden- 
tified, e.g., it is possible that the neutrino transport in 
the torus is enhanced by convective instabilities. Because 
of the quadratic dependence on the neutrino luminosi- 
ties, an increase of the neutrino fluxes would affect the 
i/p-annihilation sensitively. Alternatively, still relying on 
the simple picture described in Sects. 4.2 and 4.3, one 



might feel tempted to interpret the estimates of the an- 
nihilation energy towards the optimistic side, although 
they were derived by employing a number of very gen- 
erous and favorable assumptions. In this case some in- 
teresting implications for the hypothesis of stellar mass 
accretion disks around black holes as sources of gamma- 
ray bursts (Paczyhski 1991; Narayan et al. 1992; Woosley 
1993a, 1996) can be inferred from combining the theoret- 
ical results with information about measured burst time 
scales (Meegan et al. 1995a, Kouveliotou 1995) and en- 
ergies of cosmological bursts (e.g., Woods & Loeb 1994; 
Quashnock 1996). 

Short gamma-ray bursts have a typical duration i 7 
of several tenths of a second and a typical total en- 
ergy that is a factor of ~ 20 smaller than that of 
long bursts (Mao et al. 1994). They require the re- 
lease of neutrino energy from accretion tori with masses 
M d ~ 0.1 M [t 7 /(0.05 ... 0.2s)] -few 0.1 M and a 
beaming of the expanding fireball into a solid angle 
Sn ~ 4tt [(1 ...9) • 10 49 erg/%] [M d /O.1M ], when E 1 is 
the typical energy of a cosmological gamma-ray burst if 
the gamma-rays were radiated isotropically. <50 becomes 
smaller if the useful energy from ^-annihilation is less 
than (1 ...9) • 10 49 erg. For E 1 — 10 51 erg the focussing 
of the fireball is noticable, 8fl — 47r/10, whereas for 
Ery — 10 50 erg it is essentially absent. The required ac- 
cretion mass might suggest merging events of compact bi- 
naries as the origin of the bursts in this case. The disk 
masses and a possible focussing of the fireball expansion 
towards the observer have to be explained by theoretical 
modelling. 

For long bursts, most of which have durations t~ t 
between some 10 s and about 100 s, the accretion of 
Md — few Mq is needed and could provide an energy 
E 1 ~ (1 ... 9) ■ 10 50 [M d /M Q ] erg which is of the order of 
— 10 51 erg without significant focussing of the relativistic 
pair-photon plasma being necessary. For more energetic 
bursts some beaming of the fireball expansion would have 
to be invoked. The large accretion mass favors the "failed 
supernova" or "collapsar" scenario (Woosley 1993a). In 
this model several solar masses of material surround the 
most strongly neutrino radiating region of the accretion 



torus close to the innermost stable orbit around the black 
hole. This should lead to mixing of baryons with a signif- 
icant fraction of the pair-photon plasma and will confine 
the volume where ^-annihilation might create a relativis- 
tic fireball to a baryon-poor region along the system axis. 
The fireball will expand into a limited solid angle which 
should compensate for the reduction of the useful fraction 
of the ^-annihilation energy. 

This interpretation of the bimodal distribution of burst 
durations (Kouveliotou et al. 1993) employs two different 
kinds of astrophysical events. In contrast, Wang (1996) 
attempted to explain the bimodality by a superposition 
of two distinct time scales in the temporal structure of 
individual bursts corresponding to peak widths and sepa- 
rations between adjacent peaks. Katz & Canel (1995a, b) 
have recently suggested the association of short and long 
bursts with two different classes of models and have hy- 
pothesized that short bursts are produced by neutron star 
collisions and long bursts originate from accretion-induced 
collapse of bare degenerate white dwarfs (Dar et al. 1992). 
In both types of models the burst energy would be pro- 
vided by the annihilation of emitted vv pairs. Accretion 
induced collapse, however, was ruled out as a source of 
gamma-ray bursts situated at cosmological distances by 
Woosley & Baron (1992) on grounds of the unaccept- 
able large baryonic pollution of the surroundings of the 
collapsed star caused by a nonrelativistic neutrino-driven 
wind. The same worries also hold for collisions of neu- 
tron stars where explosions of ejected low-mass fragments 
might create an envelope of baryonic material around the 
collision site. Moreover, one has to be suspicious whether 
the neutrino emission will be luminous and long enough 
that ^-annihilation can provide an energy ^ 10 50 erg for 
a short and most likely unbeamed gamma-ray burst. 

Neutrino emission and ^-annihilation would be the 
energy source for the gamma-ray bursts also in the two 
classes of models that could lead to accretion tori around 
black holes, i.e., the merging of binary neutron stars or 
of neutron star black hole systems in case of short bursts, 
and collapsing, very massive stars which do not succeed to 
explode as type-II supernovae in case of long bursts. The 
bimodal distribution of burst durations would reflect the 
two distinct mass ranges of the accretion tori around the 
accreting stellar mass black holes, some 0.1 M or a few 
Mq, respectively. The similar peak luminosities of both 
short and long bursts (Mao et al. 1994) could be explained 
by the same underlying energy source for the gamma-ray 
bursts. The short-time variability of the gamma-ray sig- 
nal could be a consequence of jet precession (Hartmann 
& Woosley 1995) or of inhomogeneities and instabilities 
in the accretion torus that give rise to fluctuations of the 
accretion rate. And the individual characteristics of burst 
events might be associated with different masses of accre- 
tion tori and accreting black holes, different torus struc- 
tures due to different angular momentum distributions, 
and different accretion time scales because of variations 
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of the angular momentum transport, e.g., caused by mag- 
netic fields or viscosity producing dissipative processes in 
the torus. 

Movies in mpeg format of the dynamical evolution 
of all models are available in the WWW at http: 
/ / www . mpa-gar ching . mpg . de/~mor/ nsgrb . html 
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